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ABSTRACT 


This thesis consists of a three-dimensional numerical 
analysis of the Jet Aircraft Hush House located at Naval Air 
Station Jacksonville, Florida. Utilizing the PHOENICS Code 
allows ior the determination of the aerothermal 
characteristics including velocity, pressure, enthalpy, 
turbulent kinetic energy and the dissipation rate of turbulent 
kinetic energy in the facility during testing of the U.S. 
Navy's F-4 (Phantom II) J-79-GE-8 gas turbine engine with 
afterburner. How and by what method PHOENICS arrives at this 
solution is discussed. Of greatest importance is the resulting 
"behavior" of the aerothermal system. 

Problems encountered using the PHOENICS Code, resulting 
numerical solutions to the particular facility, comparison to 
actual test data and recommendations for further applications 


of the PHOENICS Code are presented. 
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ITC. INTRODUCTION 


New and overhauled aircraft gas turbine engines undergo 
stringent performance testing before returning to operational 
use. To allow for safe and effective monitoring of the engine 
and aircraft, gas turbine Test Cell and Hush House facilities 
are utilized to perform monitoring of key temperatures, 
turbine speeds and flow rates associated with the engine. The 
test facility additionally provides distortion-free air flow 
to the engine inlet (eliminating rain and wind effects), for 
lower exhaust temperatures exiting from the facility and 
reduced noise pollution to the surrounding area. Furthermore, 
these facilities allow for final engine "trim" settings to 
ensure that maximum power is available. 

Currently the test facilities designed and built for the 
Navy vary widely. The various facilities include: 


1. The Run-up pad--where a tied down aircraft conducts 
power checks in an open air environment. 


2. Sound suppressor--located at the Run-up pad to reduce 
noise. 


3. Open Test Stand (OTS)--located outdoors that supports an 
aircraft engine and has portable instrumentation 
attached. 


4. Hush House--a totally enclosed facility allowing around 
the clock testing of aircraft. 


ee fest Cell (TC)--an enclosed facility allowing testing of 
various engine sizes out of the airframe. 
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Design engineers have used experience and common sense to 
build these facilities with varying degrees of success. 
Aerothermal monitoring supplements the design of the 
SEBUGC TUES . 

With the tremendous technological advances made in recent 
years, much higher temperatures and thrusts are available 
compared to similar size and weight engines of just a few 
years ago. Considering that nearly 67% of the Navy Test Cells 
and Hush Houses are at least 25 years old it is evident that 
a new generation of standardized facilities meeting federal, 
state and local noise, pollution, security and safety 
standards is required. [Ref.1:p.3] 

The Navy Environmental Health Center cited total claims 
against the government in 1981 at $114 million due to physical 
injury, hearing, or respiratory damage to military and 
Civilian personnel. Additionally, the 1981 ruling by the 
federal court system that military gas turbine test facilities 
must comply with Section 233 of the Clean Air Act accelerated 
the need for modern test facilities meeting the above 
requirements plus those projected for the future. [Ref.1:p.3] 

The goal of the Navy is to provide for standardized 
aircraft and engine test facilities that meet both safety and 
environmental concerns while providing the latest technology 


to monitor the engines. To accomplish this, complete 


development of test facilities with reliable operation and 
Maintenance throughout the facility life cycle 1S necesSary. 
[Ref.1:p.1] 

When contemplating the inherent design criteria for the 
new facilities, problems such as: turbulence, non-concentric 
flow patterns, multiple jet engine aircraft (including varying 
temperature and thrust arrangements), supersonic velocities, 
extreme temperature and pressure gradients, and awkward 
geometries must be appraised. [Ref.5] 

Enigmas associated with the design of the structure may 
not be obvious. Complete information dealing with the 
aerothermal properties is a must in basing decisions on new 
multi-million dollar expenditures. Therefore, to minimize the 
possibility of an erroneous design, the PHOENICS Code is 
employed first on documented Test Cell and Hush House 
structures having actual field test data. After verifying the 
accuracy of PHOENICS with the test data, the Code eventually 
could be used at a Navy installation (such as the Naval Civil 
Engineering Laboratory) to design future test facilities 
providing accurate modeling data. 

In this instance, the performance of a Navy J-79-GE-8 gas 
turbine engine located at Naval Air Station Jacksonville, 
Florida Hush House facility is investigated. The data is used 


to determine the behavior of the existing structure and will 


be employed to assist engineers in the research and 


development (R&D) phase of the facility design project. 


IT. PROBLEM FORMULATION 


A. HUSH HOUSE CHARACTERIZATION 

The Hush House under investigation is an "Air Force” type 
located at Naval Air Station Jacksonville, Florida. In 1980 
the Air Force purchased 25 and later 54 demountable Hush 
Houses with Test Cell capability. This standardization 
accompanied with such a large scale purchase and the 
flexibility of aircraft or engine test capabilities allowed 
for a cost effective combined design. Because of the large 
scale purchase the cost per Hush House was a relatively low 
$2.2 million. The only major modification performed by the 
Navy to the Air Force design was inclusion of an additional 
AFFF fire fighting system. The Hush House design does not 
address air pollution standards. Some type of air pollution 
abatement system will undoubtably be required in newer 
designs. [Ref.1:p6] 

The general dimensions of the Hush House is 28.7 meters 
(94 feet) wide by 25.5 meters (83.7 feet) long with the oval 
shaped "augmenter tube” extending an additional 31.3 meters 
(103 feet) from the rear of the Hush House (See Figures 2.1 
and 2.2). The aircraft is placed so that the engine exhaust 
is approximately 4.6 meters (15 feet) from the augmenter tube 
Zaeeaice. The most common construction material for the 
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structure is steel. The design is geometrically simple 
allowing for mass producticn. The air flow in the laste 
inlets nearest the augmenter tube provides both engine intake 
and cooling air for the surrounding engine area. The final 7 
meters (23 feet) of augmenter tube is angled at approximately 
45 degrees running from the floor upward to provide for a 
smoother directional transition of gases. Acoustic damping 
materials line the augmenter tube, inlets and aircraft bay 
walls to reduce exhaust noise. 

The Hush House is used in support of the Naval Air Rework 
Facility (NARF) located at the Naval Air Station. Although the 
Hush House is presently utilized for engines fully mounted in 
an airframe, it is anticipated that engine "Test Cell” testing 
(where an engine is mounted on atest stand) will be approved 
in the future. The facility is normally in operation from 0700 
to 2400 Monday through Friday with an average of 1 to 5 
aircraft being tested per day. One of the most common aircraft 
tested is the F/A-18 Hornet powered by two General Electric 
F-404 gas turbine engines. Typically, after completicnmion 
overhaul of the aircraft and engines, the aircraft will be 
towed to the facility and backed into the Hush House. The hugh 
aircraft inlet doors to the facility will then be closequaee 
provide an enclosed area free of foreign matter, outside wind, 


rain, and other environmental and personnel concerns. 


With the aircraft fully enclosed and securely tied down, 
a normal crew of three mechanics will perform final engine 
"trim". One person is located in the Hush House Control Booth 
to ensure manning of the fire fighting systems. The second 
will physically adjust the fuel control mechanism on the 
engine via radio signals received from the cockpit where the 
third crewman is. The engines are controlled during testing 
by the third crewman in the aircraft who is additionally 
responsible for monitoring the various fuel, oil, % power, 
hydraulics, and other gauges. Testing normally lasts for a 
period of two hours. Only one engine is tested at a time with 
the other normally at idle. Air intake and cooling is provided 
by five inlet areas located on either side of the Hush House 
(See Figures 3.2 and 3.3). Wire screening, air diffuser and 
noise abatement equipment 1s present at the inlets. Installed 
fire fighting and wash-down equipment is present in the 
a@ircraft bay area with the air inlets capable of being 


automatically closed. 


B. ISSUES CONCERNING THE HUSH HOUSE 

The Hush House is currently performing as designed with 
fee Delow exceptions. Investigation indicates some cracking 
of the welds which fuse together the stainless’ steel 
“corrugated" sheet metal overlaying the sound absorption 


feet ial located in the augmenter tube. Also, actuation of 


heat sensors in the overhead of the aircraft bay sometimes 
occurs when the engines are idled for a period of time. 
Apparently some of the exhaust gases travel upward when low 
mass flow rates are present instead of going down the 31 plus 
meters (103 feet) of augmenter tube. 

In general, personnel at the Naval Air Station were 
pleased with the overall operation of the Hush House. Although 
vibration of the entire physical structure has been noted in 
this facility as in other designs, it does not appear to cause 
any structural problems such as the requirement for inspe@uumen 
and retightening of bolts as occurs in other Navy Hush House 
designs. 

Physical test data is available for an identical Air Force 
facility for the F-4 Phantom II and is included in partcues 
this report for comparison with results from the PHOENICS 
Code. The test data provides an important comparison for this 
thesis to verify the computer model. [Ref.2] 

A physical inspection of an identical Hush House to Naval 
Air Station Jacksonville's was under-taken in June 1988 at 
Dobbins Air Force Base in Atlanta Georgia. Discussions with 
the various Air Force personnel corroborated the previously 
mentioned problems with the exception of the heat sensor 


activation problem. The Dobbins Air Force Base Hush House was 


Memirzead for both Test Cell and Hush House operations 


Mmevwolving the F-15 (EAGLE) fighter . 
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Figure 2.1 Actual Hush House (Top View Looking Down) 
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Figure 2.2 Actual Hush House (Side View) 
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IIT. THE PHOENICS CODE 


A. AN EXPLANATION 

PHOENICS, which stands for "Parabolic Hyperbolic Or 
Elliptic Numerical Integration Code Series", is a computer 
code capable of simulating fluid flow (including changes of 
state), heat transfer and chemical reactions. The code is 
copyrighted by CHAM (Concentration Heat and Momentum) 
Limited which was founded by D. Brian Spalding in 2997. 
(Ret. pata) 

As shown in Figure 3.1, the PHOENICS Code is structured 
into three major components consisting of “BARToee 
"SATELLITE", and "GROUND STATION". The "EARTH" program may be 
considered as a "black box" which contains the computational 
procedures for all applications. The "SATELLITE" program 
consists of several data-input methods which define the 
boundaries and scope of the problem. Finally, the "GROUND 
STATION" program haS various subroutines containing built- 
in physical models which may or may not be selected depending 
on the program application. 


By default "EARTH" can solve for up to 25 dependent 


variables. In the Hush HousSe application these include 
pressure, all three dimensional velocity components, 
enthalpy (temperature), turbulent kinetic energy, and the 
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turbulent kinetic energy dissipation rate. The dependent 


variables can be enlarged if desired. 


| 
PHOENICS input library 


at P, 


Satellite 


ed 





Figure 3.1 Satellite, Earth and Ground 


The PHOENICS Code is exceptionally flexible allowing for 
varying sizes and scales of structures from "micro" to 
"Macro". A "“final" solution is highly dependent on the 
accuracy required by the user as well as the grid sizing used. 
The higher the accuracy Biome manbes of "cells" desired, the 


greater the computational effort required for a solution. In 
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this investigation, a "grid independent solution", where a 
sequence of "finer" grids would be used, was not possible due 
to time limitations. 

In order to solve a given problem the user must create a 
"Ol" file. The "Q1" file describes all aspects of the system 
being investigated. This file is partitioned into several 
groups allowing the user to specify various aspects of the 
physical system. Examples are structure size and shape. 
Additionally, other groups are capable of providing the 
necessary grid Spacing, boundary conditions, engine mass 
flow rate, engine exhaust temperature and printing options. 
The "Ol" file is read by the "SATELLITE" program which is 
recognized and processed by "EARTH". The "EARTH" program 
solves discretized conservation equations. A typical 
conservation equation (before discretization) which "EARTH" 


solves can be expressed in the form: 


- (rp@) + -divirpve? -— ris graded) s—s.. 
o£ 
where: 
r = phase volume fraction 
@ = dependent variable 
Tes = exchange coefficient (laminar or turbulent) 
S+ = source or sink term 
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density 


O 
m 


< 
UT 


velocity vector 


Here @ is the dependent variable currently in question. For 
example, using the conservation of energy equation, @ would 
be enthalpy and [fs is the diffusion constant. [(Ref.3:p.2.4] 
The values computed and displayed for each cell are the 
solutions of the algebraic equations relevant to the laws of 
physics for that cell. The laws of physics are those of 
conservation of mass, energy and momentum and any other 
dependent variable. The general form of the algebraic 
linearized equations resulting from a discretization of the 


previous field equation is: 


ace + awhw + Ann + Aasds + Aanpda + ardr + ardGr + S 


@e + aw +t a@xn +t a@s t+ Aan + At + A@1 + Ap 


where ¢@ 1S the value of the dependent variable at a 
particular cell and subscripts E, W, N, S, H and L are 
neighboring locations. The values of $ at these locations are 
known from a previous iteration. Quantities S and ap express 
the influence of a source of the entity ¢. The a's are 


influence coefficients, temporarily treated as constants, with 


des 


the above subscripts expressing interaction between neighbor 
cells. [Ref.32p.2.43 

Solvang the multitude of algebraic equations is 
complicated and time consuming. An iterative manner is used 
until the imbalance between the left and right sides of all 
the above equations are negligible. Several schemes can be 
used to perform the "Sweeps" (where a sweep solves for the 
dependent variables for each cell in the grid one time) 
necessary for the iteration process. The results will be 
identical although the time to complete the task will vary 
with the iteration scheme. [({Ref.3:p.2.5] 

There are several ways to determine whether the problem 
has "converged" to a solution. Results are monitored to 
determine whether values of dependent variables are remaining 
constant or fluctuating. The "error" (residuals) is tabulated 
for each sweep at the end of the computer printout. By 
observing "residuals" from each sweep a pattern of decreasing 
("converging") values should be present. If the "residuals" 
are increasing ("diverging") adjustments of several constants 
which in-turn apply a dampening effect upon the variables in 
question may be required. The printout of successive residual 
values demonstrate convergence by a reduction in the 
residuals. That is, a visually downward sloping line on the 


accompanying graph Of residuals shows convergence. 
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Conservation balances on each dependent variable can be 
verified using "net source" values which are tabulated at the 
end of each run. The mass flow rate (mass flow rate = density 
* area * velocity) 1s additionally calculated, in this 
instance, for the five inlets and one outlet confirming that 
no hidden, missing, or improper boundaries are present. The 
problem has "converged" when an acceptable level of error 
remains. 

Successful “convergence" of the PHOENICS program is 
dependent upon experience and an understanding of the results. 
If the residuals of a dependent variable show divergence, the 
user can manipulate the rate of change of values of that 
variable by using the RELAX command to return the variable 
back to the convergence path. Relaxation is placed on a 
variable to control the variables maximum change fron each 
sweep. This minimizes the possibility of "over-correction" 
or “blowing up" of the values. Smaller values of relaxation 
result in smaller value differences per iteration. However, 
excesSively small relaxation values require additional 
computer time. A balance gained by user experience enables 
optimal convergence with minimal computer time. 

When multiple runs are required for large programs output 
data 1s stored ina "“DFO9" file which allows the next run to 


resume the computations (i.e., a restart) where they last 
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ended. The complexity of the problem and the required 
accuracy determine the quantity of computer runs needed. The 
number of runs can vary from the tens to the thousands and 


higher. 


B. ASSUMPTIONS IN THE MODE, 

There are several approximations that deserve mention in 
the modeling of the Naval Air Station, Jacksonville, Florida 
Hush House. Because the geometry of the model uses cartesian 
coordinates, the actual Hush House's rounded and curved 
surfaces, including the augmenter tube, Hush House overhead, 
engine and exhaust gas exit ramp, were approximated as flat 
stepped surfaces. 

The inlet "muffler" arrangement has been eliminated for 
Simplification. The corrugated metal found througnhcutweene 
structure to cover the sound absorption material was 
considered to have a friction surface. The model, therefore, 
included friction for all walls, flocrs7 ceca. 

The actual test data used for comparison between the 
model and facility contain several notable differences. The 
model uses the symmetry of the centerline of the Hush House 
to reduce the number of cells by half. The engine is located 
on the centerline in an ideal cradle position. Because of 
the symmetry only half of the engine is observed as if the 


Hush House was being used as a test cell without the aircraft 
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Bresenc. Since the correct engine exhaust area is used the 
units of mass flow rate and temperature are correct for the 
model configuration. The engine is configured for afterburner 
operation. The physical aircraft is not present. 

Actual test data used for comparison consists of the 
entire aircraft being placed in the facility. The engine 
configuration has both engines running with one at 
afterburner while the other engine is at "engine idle". The 
obvious result of the engine arrangement 1s a non-symmetrical 


engine location in comparison to the symmetric model. 


Cc. MODELING THE PROBLEM 

Data, including blueprints of the Hush HouSe and aircraft 
engine information, was gathered to accurately depict the 
facility configuration. The information obtained was in 
English units and was converted to the International Metric 
Byectem Or "S.1I.”" using standard conversion factors. 

With the appropriate information gathered it was time to 
"fit" the Hush House to the cartesian coordinate system of X, 
Y and Z coordinates using PHOENICS version 1.3. The cartesian 
coordinate system was choosen because of its simplicity. Since 
the actual structure does have rounded surfaces some 
approximations were made specifically to the augmenter tube 
and overhead of the inlet bay. Additionally, the aircraft was 


eliminated to reduce the complexity of the problem allowing 
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the structure to serve in its "Test Cell” Configurar ues 
Drawings are found in Figures 3.2 and 3.3 of this se@euen 
illustrating the key items found in the structure. The 
drawings view through the structure to note the key items, 
such as: air inlets, augmenter tube and engine locations. 
Additionally, the structure walls were considered to be of 
infinitesimal thickness to further simplify the grid 
structure. The engine was simplified from a cylindrical shape 
to one of a rectangle of equal intake and exhaust area as the 
original cylinder. Other small items, such as the engine 
stand, were eliminated. 

To take advantage of the symmetry of the Hush House and 
thereby reduce computer CPU time, only half of the Hush House, 
engine and augmenter tube were modelled in the "901" file. This 
was easily accomplished by dividing the structure at the 
centerline. This also appreciably reduces the size of the 
@Oo tat i ler 

Figures 3.2 and 3.3 graphically display the X-Z anG@iiene 
Y-Z planes of the structure. Intuition, trial-and-error yea 
available computer time dictated the grid spacing used. 
Obviously, volume areas of interest such as the engine exhaust 
and augmenter required smaller cell sizZe. 

The "Ol" file describes the structure using several types 


of statements. The "CONPOR" statement allows the user to 
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define WallovowmernNeadcusanmces the stioor- The “PATCH” and 
EPeCOVAI” statements work together to indicate items such as 
an engine exhaust or atmospheric boundaries. Dependent 
variable values for say the atmospheric boundaries can be 
"fixed" or allowed to "float" depending on the problem 
definition. When the walls, inlet bay overhead, augmenter 
tube, engine and floor dimensions have been measured, a grid 
for the system can be constructed, based on locations and 
anticipated responses over different "areas" of the system. 
Areas of anticipated large response gradients would have 
smaller cell volumes than low activity areas. 

Netemthatmin figure 3.2 the Z-axis consists of 40 cells 
consecutively numbered from 1 to 40 going left to right. The 
X-axis consists of 10 cells consecutively numbered 1 through 
Meegoing from lower to upper. Figure 3.3 is similarly numbered 
Biere the Y-axis rune from 1 to 12 going from lower to upper 
and the Z-anis is the same as previously described. 

iemsacecmitdon to Geccrabing the geometric structure and 
Prescribing the "PELAXATION" devices, the properties of the 
fluid medium must be specified. The following boundary 
conditions were used: 

1. ambient temperature = 295 K 
2. ambient pressure = 101325 Pa 


3. gas constant = 286 Joules/(kg —- R) 


OBA: 


4. Cp = 1004 Watt-sec/(kg ser) 


5. mass flow rate (engine) = 197.72 kg/(sec-sq meters) 
6. temperature (engine out) = 1944 K 
Additionally, "friction" was applied to all surfaces of 


the Hush House. Appendix A contains the final form of the 
Q1 file for reference. An in-depth explanation of all the "Qi" 
program parameters can be found in the “Beginners Guide and 
Reference Manual" (See Reference 4). 

Several methods and theories have been previously used to 
attempt to reach a "“converged" solution in other test 
facilities. The method chosen for this project was to 
initially minimize the “excitation” of the model thereby 
reducing the "jump" required in solving the momemau7e 
enthalpy, turbulent kinetic eneray dissipation rate. turbulent 


kinetic enerscy and mass equations. To accomplish this, during 


Successive runs (where a "run" consists of several “Ssweeps"), 
small incremental increases to the source were applica es 
jet engine exhaust temperature and mass flow rate were slowly 
increased, from run to run, while carefully monitoringeeme 
residuals as well as the mass flow rate in and out of the 
Structure. the eurbulenee model, "KEMODL", was used 
throughout rather than considering the flow first as laminar 
and then adding turbulence. This method appears to have met 


with considerable success. The table that follows illustrates 


the iterative process used throughout the computer modeling: 


vy) ery 
on 


Taplemeomr ILTERATIVE PROCESS 


! Temperature | Flow Rate | Sweeps | 
| 

i 400 f 10 f 600 | 
f 2101.6, A 25 , 600 | 
; 750 i 25 [ 600 ] 
f 750) f 50 ; 600 ! 
f HOO | fies) i 2500 
i 1250 | eS A E200 i 
f 250 | 100 i L200 fH 
| JL 4/518) | £25 | 1200 | 
| 1500 | 25 | ZO | 
| fs ie) ; 2 © | 2000 | 
| aes) | AS: | 2000 | 
| 1944 | SiO | 2000 | 
1944 | nee) | 5000 | 
| 1944 | m7. 7 2 | 5500 
ee ee ee 

| Total Sweeps = 26,200 | 


EE ————EE—————E—— ee | 


This amounts to approximately 147 hours of computer CPU time 
on the IBM 3033 mainframe computer installed at NPS. 

Using the MVS “batch" system on the IBM 3033 mainframe 
proved to be time consuming. The PHOENICS Code is designed to 
work both interactively and in batch mode. Due to the high 
computer usage at NPS the "interactive" mode was not 
available. Limited to the batch mode, careful planning for 
evening or night runs was required to maximize CPU time. The 
number of "Degrees of Freedom" (DOF) that must be solved for 
gives an idea of the problem complexity. DOF = NX * NY * NZ 


Seene 7 Gependent variables. Solving this for the grid size 
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of 10 * 12 * 41 yields nearly 35,000 values calleulatecaduiaas 
each sweep. Typically, 150 to 200 sweeps could be accomplished 
during the day using four 15 minute (class "G") runs GEM 
sweeps each. After normal working hours three one hour (class 
"J") runs totalling 500 to 600 sweeps were accomplished. Data 
stored on the all important DFO9 file was continuously 
"backed-up" after a successful run to prevent loss of data if 
an improper subsequent run resulted in erroneous data. The 
backed-up data could then be restored to the DFO9 file to 
enable continuation of the runs without starting at the very 


beginning. 
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Figure 3.2 Hush House Model (X-Z plane) 
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Figure 3.3 Hush House Model (Y-Z plane) 
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IV. RESULTS 


A. THE FINAL SOLUTION 

After obtaining the "solution" for the problem the 
question still remaining is just how "good" are the computer 
model results? Some of the "tools" available for use, such as 
checks on reSidual values and mass flow rate calculations, 
were discussed in Chapter III. The checks indicate that a 
converged solution is obtained, but, still do not prove the 
computer model data validity with the actual facility data. 

Using actual (experimental field measurements) results 
from Reference 2 allows for an accurate evaluation of the 
computer model. Appendices "C" and "D" offer some experimental 
data for velocities and temperatures, thereby indicating mass 
flow rate and augmentation ratio. Additionally, figures 4.1, 
4.2 and 4.3 found in this section give the numerical 
(computer) model results of differential pressure, temperature 
and axial velocity profiles encountered at the IX=1 slab. 
Graphics generated pictures (generated by "PHOTON" which 
converts the computer generated results to graphical format) 
allow for the interpretation and digestion of the nearly 300 
pages of numerical output of a PHOENICS run. A condensed 
output is found in Appendix "B". In this case, PHOTON graphics 
Specifically provides graphical information for a better 
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interpretation of velocity fields and both pressure and 
temperature profiles and contours enabling a better 
understanding of computer generated data. Use of all these 
additional "tools" allow for a determination as to whether the 
computer results are valid. 

Observing figures 4.4 and 4.5 (top looking down) show the 
flow pattern of the air entering the five air inlets into the 
Hush House at IY=4 (engine level). A "typical" low velocity 
air flow region (Figure 4.4), illustrated by the minuscule 
length of the velocity arrows (velocities of less than 2 
meters/sec or 6.5 feet/sec), is noted between the first two 
and final three inlets. Figures 4.6, 4.7 and 4.8 illustrate 
the side view velocity flow patterns for the overall facility, 
engine area, and exhaust to atmosphere respectively. Note the 
ever reducing velocity magnitudes which vary from 818 meters 
per second (2684 feet/sec) at the engine exhaust to less than 
100 meters per second (328 feet/sec) as the gases travel 
further away from the engine and down the augmenter tube until 
reaching the exhaust. 

Figures 4.4, 4.5 and 4.9 (top looking down) show the air 
and gas flow patterns entering and exiting the engine region 
and augmenter tube at engine level. Figure 4.10 indicates the 
flow patterns at the top of the augmenter tube (IY=8) for the 


entire Hush House. Of interest is the counter-flowing currents 
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at the top of the augmenter tube starting at the engine 
exhaust and running to nearly half way down the augmenter tube 
ecself. 

The engine occupies cells defined by IX=1, IY=4 to 5 and 
IZ=12 to 14. As anticipated, figures 4.4, 4.9 and 4.11 show 
that at engine level the air in the Hush House moves toward 
the engine inlet from the surrounding area while air also is 
entrained at the exhaust of the engine. Figure 4.1 illustrates 
this entrained "effect" noted by the high differential 
pressure of the exhaust gases immediately exiting the engine 
at IX=1 and IY=4. Less than a meter after exiting the engine 
exhaust, the exhaust gases and inlet air have begun mixing 
which results in an area of drastic temperature reduction for 
the exhaust gases in the Hush House and immediate augmenter 
area. Gradually, as the gases pass further downstream, 
pressure begins rising which reduces the velocity of the gases 
in the augmenter tube until the gases exit to the atmosphere. 
At the inlet area to the engine (at Z=18 meters on figure 4.1) 
a low pressure area 1S noted as one would also anticipate as 
the air is drawn into the engine inlet. 

The exhaust gases exit the engine at 1944 Kelvin (3040 F). 
The inlet air and exhaust gases mix due to the entrainment 
"effect" causing a distinct drop (See figure 4.2) in 


temperature as the gases travel further down the augmenter 


ag 


tube. Temperature contours at the IX=1 slab shown in figures 
4.12 and 4.13 portray the temperature patterns out of the 
engine progressing down and out of the augmenter tube. 
Analogous temperature contours at the IX=3 slab are presented 
in figures 4.14 and 4.15. Figures 4.16, 4.17 and 4.18 empea 
the temperature contours at the IY=5 (engine) level while 
figure 4.19 shows the temperatures at the top of the augmenter 
tube at Iy=8. The temperature contours indicate a 
progressively widening area of decreasing exhaust gas 
temperatures where the vast majority of the temperature drop 
occurs in the first five meters after the engine exhaust (the 
temperature differential goes from 1944 K to 800 K). 
Differential pressure contours of the engine exhaust at 
IX=1 are shown in figures 4.20 (IZ= 14 to 24) and 4.21 (IZ= 
25 to 40). As indicated in the pressure profile, the exhaust 
jumps to nearly 6000 pascals above ambient upon exiting the 
engine, but soon after drops into a low pressure region until 
reaching ambient pressure half way down the augmenter tube at 
IZ=24. Gradually, the pressure builds until the gases exhaust 
to the atmosphere. Similar contours for the IX=3 slab are 
included in figures 4.22 and 4.23. Pressure contours foreene 
X-Z planes at IY=5 (engine level) and IY=8 (top of the 
augmenter tube) are contained in figures 4.24, 4.25 and 4.26. 


Of noted interest is how the differential pressures remain 


30 


fairly uniform at a given IZ location as the gases travel down 
the augmenter tube with the exception of the region initially 


downstream from the engine exhaust. 


B. ACTUAL VERSUS MODEL COMPARTSONS 

Generally, the results from the computer model closely 
approximated the actual. The below "actual" data, and that of 
Appendix "C" and "D" for the computer model are indicated in 
Been "English" and "S.I." units below. A comparison of the 
PHOENICS output visually shown in Appendices "C" and "D" 
illustrates the dissimilarity between the actual and computer 
model. As indicated by the data, the differences in the 
computer model (235.5 ft/sec or 71.8 meters/sec) and actual 
(227.2 ft/sec or 69.3 meters/sec) averaged velocities found 
in the augmenter tube at slab IZ=27 differed by less than 4 
percent. Averaged temperatures in the same location differed 
between 401.6 degrees Fahrenheit (47& kK) for the computer 
model and 454 degrees Fahrenheit (507 K) actual with a 
resultant difference of 6 percent. Reference 2 does not 
indicate whether the effects of radiant energy from the engine 
exhaust are considered for the actual temperature readings 
taken in the augmenter tube. If these effects were not 
considered the 6 percent difference would be reduced. The Hush 
House mass flow rate for the computer model (2373 lbm/sec or 


1076 kg/sec) versus the actual (2067 lbm/sec or 938 kg/sec) 


om 


differed by 14.8 percent. The computer model mass flow rate 
value was actually doubled (hence the 2373 lbm/sec or 1076 
kg/sec value) to "model" the entire engine and Hush House. 


The Augmentation Ratio is defined as follows: 


Mine: GS Mengine 
A.R. = ---------------- 
Miniet 
where: 
Miniet = inlet mass flow rate 
Mengine = engine mass flow rate 
The inlet mass flow rate = air density * five inlet 


opening areas * averaged inlet velocity. The engine mass flow 
rate is given at 197.72 kg/sec/sq meters * 0.195 sq meters 
engine exhaust area for half of the engine. 

Augmentation ratios range from 5 to 15 depending Ones 
Gesign of the Hush House or Test Cell where the higher the 
ratio typically indicates a more efficient design. The 
computer model augmentation ratio was found to be 10.65 while 
the actual was 8.3. The resultant difference is therefore 
nearly 28 percent. Considering the limited data taken at this 
actual power level (See Appendix C) it is quite conceivable 
that the actual augmentation ratio could be greater. 

For the actual results the velocity, pressure and 


temperature measuring devices were photographed because of the 


ie 


large fluctuations of the gauges. Additionally, velocity data 
at the final three inlets could not accurately be taken making 


it difficult to accurately determine the augmentation ratio 


and the true velocities at the inlets. 
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4.3 Axial Velocity Profile at Volume Slab IX 
for IyY=4 
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Figure 4.4 Flow Vector Plot of Hush House Inlets and 
Engine Area (X-Z Plane at IY=4) 
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Flow Vector Plot in X-Z Plane at IY=4 


(Engine Exhaust and Augmenter Tube) 


Figure 4.5 
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Flow Vector Plot in Y-Z Plane at IxX=1 


Figure 4.6 


(Entire Grid) 
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Flow Vector Plot in Y-Z Plane at IxX=1 
(Engine Exhaust Area) 


Figure 4.7 
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in Y-Z Plane at IX=1 
(Exhaust to Atmosphere) 


Figure 4.8 Flow Vector Plot 
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Flow Vector Plot in X-Z Plane at IlY=4 


Figure 4.9 


(End of Augmenter Tube) 
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Figure 4.10 Flow Vector Plot in X-Z Plane 
(Mirrored Image of Entire Grid) 
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Figure 4.11 Flow Vector Plot of Augmenter Tube Area 
for X-Z Plane at IY=4 
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Temperature Contour in Y-Z Plane at IxX=1 


(Engine and Augmenter Tube) 


Figure 4.12 
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Figure 4.14 
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Figure 4.16 


(Mirrored Image of Engine Exhaust) 
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Temperature Contour in X-Z Plane at IY=5 


Figure 4.17 


(Engine Exhaust) 
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Figure 4.18 
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Figure 4.19 Temperature Contour in X-Z Plane at IyY=8 
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V. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

Using the Parabolic, Hyperbolic, or Elliptic Numerical 
Integration Code Series (PHOENICS) a computer model of the 
Naval Air Station Jacksonville Florida's Hush House facility 
with a J-79-GE-8 gas turbine engine was simulated. The 
performance of this computer model was compared to actual 
performance data taken with an F-4 Phantom II undergoing 
testing at an identical facility located at Kelly AFB. An on- 
site inspection was conducted at Dobbins AFB in Atlanta 
Georgia to provide "hands-on" experience with the facility. 

The computer analysis consisted of slowly building up 
excitation in the model (using the KEMODL turbulence model) 
until full power was obtained. Computer generated (tabulated 
and graphical) results were compared with the actual test 
data. Results of the computer model, in comparison to the 
actual data, were in good agreement. Flow directions, 
velocities and temperatures were consistent with actual 


results. 


B. RECOMMENDATIONS 
A comparison of actual to computer model data for this 


thesis indicates that PHOENICS is a useful design and 
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troubleshooting tool allowing the user to better understand 
the physics of the model. PHOENICS at NPS 1s now just coming 
into potent use with a variety of research projects for 
aeronautical and mechanical engineers. 

To ensure proper and full use of this code several major 
"glitches" must be overcome. Time constraints placed on the 
use of the IBM 3033 main-frame computer limit the number of 
"sweeps" obtainable by the student. The batch system does not 
lend itself to correcting possible errors that may arise. An 
interactive system, such as the "VAX" stand alone systems 
being installed in Mechanical Engineering, is an ideal answer. 

PHOTON was not available with the PHOENICS installation 
on the IBM 3033 main-frame computer. A complicated conversion 
process (See Appendix "E") was undertaken to convert IBM to 
VAX "format" to allow the use of PHOTON using data from the 
Main-frame. This conversion process which involved a swap of 
information from IBM to DEC to VAX hardware took nearly 2 
weeks to accomplish. In September 1988 this transfer procedure 
will no longer be possible at NPS with the elimination of 
certain hardware at the Computer Science Department although 
data transfer will still be possible to other facilities such 
as CHAM. 

The use of PHOTON early in the computer "runs" would have 


uncovered several missing "CONPOR" statements and therefore 
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“verified” the Structural integrity of the facility 
Additionally, the amount of time required to "wade" through 
the 300 pages of printout could be reduced to several video 
screens with a much reduced volume of printout. 

The computer model used in this thesis was formulated 
using cartesian coordinates. Several parts of the structure, 
including the augmenter tube and engine, were curve shaped. 
Geometric approximations were made to substitute for the curve 
shaped areas. A more computationally time consuming model with 
increased numbers of cells and/or the use of body-fitted- 
coordinates (BFC) would enhance the computer model. 

As noted in Figure 4.2, over fifty percent of the 
reduction in engine exhaust gas temperatures occurred in the 
first 20 percent of the augmenter tube. Could an improved 
design include shortening the augmenter tube? This type of 
"what if" question can be answered using PHOENICS. With the 
variable data file ("DF09") already in place the grid pattern 
can be modified to shorten the augmenter tube. Minimal 
computations would be required to gather definitive results. 

Parameter studies, including shortening the length of the 
augmenter tube (to determine the effect of augmenter tube 
length on performance), movement of the engine location in the 
Hush House, changing the inlet locations and sizes, and the 


use of different turbulence models could lead to improvements 


on the present Hush House design. In fact, the ability to 
change parameters and observe their effect on system 
performance is the prime attribute of numerical modeling. 
Finally, a discussion of turbulence is in order. PHOENICS 
offers four turbulence models which can be readily selected. 
The four models are: Algebraic, Mixing Length, K-Mixing Length 
or K-Epsilon. The user can also formulate a model of his own 
choosing which can be incorporated into the code via the 
"GROUND" program. Two questions arise. First, which turbulence 
model is appropriate for the problem under investigation? 
Secondly, how strong is the coupling between other system 
response variables and turbulence variables? The second 
question might be answered by running the PHOENICS Code on the 
problem using several different turbulence models. 
Unfortunately there was not adequate time to perform this in 
the investigation. The first question can not be resolved 


until a better understanding of turbulence is obtained. 


Cc. CONCLUSIONS 

The computer model results generated by this thesis 
closely parallel the actual test data. The results therefore 
indicate that the PHOENICS code can be successfully utilized 
to model future Test Cell and Hush House applications as well 
as “trouble-shoot”" suspect areas of interest. Also, areas that 


may be difficult to properly monitor with test equipment could 
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be modelled instead. With the use of additional “cells" and/or 
BFC a more refined model could more closely duplicate the test 
facility providing more accurate results. Follow-up parametric 
studies to improve and verify the design and data seem in 
order. The results indicate that the Hush House is performing 


as intended. 
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APPENDIX A 
THE PHOENTICS ©1 FILE 


TALK=F ,RUN( 1,1) 
FROG III RO IER RIITIGE IIIB RE UIE mar at an 
# NAVAL AIR STATION JACKSONVILLE FLORIDA HUSH HOUSE TEST * 
* FACILITY. THESIS PROJECT OF LT ERIC A. NICOLAUS USCG~ «# 
* WITH THESIS INSTRUCTOR PROF DAVID SALINAS. * 
JSR Sr Si BE RET OE IIIS SHE A el 


URIS ORBIT HRI EIHIRI HIRE RIMM OI IBHIHRBRIRFIOH NI AHI IR HHHEE ne 
LIST OF VARIABLES 


MJET: MASS FLOW RATE OF GAS (KG/SECOND/SQ METER) 

TJET: JET EXHAUST GAS TEMPERATURE (DEGREES KELVIN) 

PJET: EXHAUST PRESSURE OF JET ENGINE (NEWTONS/SQ METER )* 
RHOJET: DENSITY OF THE GAS (KG/CUBIC METER) 

HJET: ENTHALPY (JOULES/KG OR NEWTON*METERS/KG) 
KEJET: KINETIC ENERGY (JOULES OR NEWTON-METERS) 
GASCON: GAS CONSTANT (JOULES/(KG-DEG KELVIN)) 
TAMB: AMBIENT AIR TEMPERATURE (OEGREES KELVIN) 
HAMB: AMBIENT ENTHALPY (JOULES/KG) 

KEINJ: TURBULENT KINETIC ENERGY (JOULES) OF JET 
CSUBP: SP HT 2 CONST PRE (JOULES/(KG-DEG KELVIN)) 
EPINJ: DISSIPATION RATE OF TURBULENT KE (JOULES) OF JET 
KEINA: TURBULENT KINETSC ENERGY (JOULES) AMBIENT 
EPINA: DISSIPATION RATE OF TURBULENT KE (JOULES) AMB 
Pl: PRESSURE (NENTONS/SQUARE METER) 

Hl: ENTHALPY (JOULES/KG)} 

Ul: VELOCITY IN X-DIRECTION (METERS/SEC ) 

V1: VELOCITY IN Y-DIRECTION (METERS/SEC ) 

W1: VELOCITY IN Z-DIRECTION (METERS/SEC ) 

TMP1: TEMP OF THE FIRST PHASE (DEGREES KELVIN) 

RHO1: DENSITY OF THE FIRST PHASE (KG/CUBIC METER) 
PRESSO: REFERENCE ATMOSPHERIC PRESSURE (NEWTONS/SQ mir )* 
TMP1A: GROUND REF TEMP USED IN GRND2 (DEG KELVIN) % 
TNHP1B: I/7CSUBP USED IN GRND2 EQU (KG-DEG KELVIN/JOULES )* 
ORHIDP: DELTA RHO/DELTA PRESSURE W/RESPECT TO H1 # 
FACT: CORRECTION OR "FUDGE" FACTOR % 
ENUL: LAMINAR KINEMATIC VISCOSITY (METERS SQ/SEC) * 
ENUT: TURBULENT KINEMATIC VISCOSITY (METERS SQ/SEC) % 
96 EF III IEE IR OE III EERIE III OTE 


* 
a 
WJET: VELOCITY OF ENGINE EXHAUST (METERS/SEC ) # 
*% 
% 


xx xk kK KK 


xee eK KK KK KKK KKK KKK KKK KKK KK KK KK SK 
xk «x kKKK KK SY 


GRCUP 1. RUN TITLE AND OTHER PRELIMINARIES 
TEXT(NAS JACKSONVILLE HUSH CHAMBER ) 
REAL(WHJET ,»MJET,TJET ,PJET »,RHOJET »,HJET ,KEJET ,GASCON,CSUBP , TAMB ) 
REAL( HAMS ,KEINJS, EPINJ,KEINA,EPINA ) 


GROUP 2. TRANSIENCE; TIME-STEP SPECIFICATION (DEFAULT=STEADY ) 


GROUP 3. X-DIRECTION GRID SPECIFICATION 
NX=10 
XFRAC(1)=5-1 35 XFRAC(2)20.31 
XFRAC(3)=G 5 XFRAC(4}=0.5975 
XFRAC(5)=35 5 XFRAC(6)=2.0 
XFRAC(7)22 3 XFRAC(8)=1.9 
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GROUP 4. Y-DIRECTION GRID SPECIFICATICN 
NY =12 
YFRAC(13=-3 5 YFRAC(2)=0.562 
YFRAC(3 )=2 § YFRAC(G )=0.314 
YFRAC(5 )=3 § YFRAC(6)=0.562 
YFRAC(7}3=4 s YFRAC(8&)=0.75 


GROUP 5. 2Z-DIRECTION GRID SPECIFICATION 


NZ=40 

ZFRAC(1 J=-1 ZFRAC( 2 }=1.0 

ZFRAC(3)=1 ZFRAC(4)=1.7 

ZFRAC(53=1 ZFRAC(6 }=2.0 

ZFRAC(7}=1 ZFRAC(83=1.7 

ZFR&C(93=2 ZFRAC(10)=3.0 


ZFRAC(11)=1 
ZFRAC(13 )=1 
ZFRAC(15 )=1 
ZFRAC( 17 3=2 
ZFRAC(19)=1 
ZFRAC(21)=2 
ZFRAC( 23 )53 
ZFRAC(253=2 
ZFRAC(27)=7 
ZFRAC(29)=1 
ZFRAC(313=4 
ZFRAC(33354¢ 
ZFRAC!{ 35 )=4 
ZFRAC(37)=1 


ZFRAC(12)=1.5 
ZFRAC( 14 )=1.0 
ZFRAC(16 3=1.7 
ZFRAC(183=1.0 
ZFRAC(20)=1.7 
ZFRAC( 22 3=1.0 
ZFRAC( 24 3}=0.566666 
ZFRAC(26)=0.35 
ZFRAC( 28 3)=3.0 
ZFRAC(30 )=2.0 
2FRAC(323=0.75 
ZFRAC( 34 )=0.575 
ZFRAC(36)=0.5 
ZFRAC(38)=3.0 


wo ww wee wwe wee ee 


GROUP 6. BODY-FITTED COORDINATES OR GRID DISTORTION 


GROUP 7. VARIABLES STORED, SOLVED & NAMED 
SOLVE(P1,H1,U1,V1,W1) 
SOLUTNIP1,Y,Y,Y,N»N»N) 

STORE( TMP1,RHOI ) 


GROUP 8. TERMS (IN DIFFERENTIAL EQUATIONS) & DEVICES 


GROUP 9. PROPERTIES OF THE MEDIUM 
TAMB=295. 
TJET=1944¢. 
PJET=0. 
MJET=197.72 
PRESSO=101325. 
RHO1A=0. 
GASCON=286. 
RHO1B8=1./GASCON 
RHOJET=( (PRESSO+*PJET }®#RHO1B )/TJET 
WJET=MJET/RHOJET 
TMP 1A=TINY 

*e% TMP1B = 1. DIVIDED BY C-SUB-P3} USED IN TMP1=GRNDS5 

CSUSP=1004. 
HAMB=TAMB*CSUBP 
TMP 1E=1./CSUBP 
HJET=CSUBPRTJET 
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we NEXT 2 STATEMENTS INVOKE COMPRESSIBLE FLOW IN GREX1 
TMP 1=GRND2 
RHO1=GRNDS 
DRH10P=GRNOS 
ENUL=1.E-05 
ENUT=200. *ENUL 
TURMOOI( KEMOOL ) 
STORE(VIST) 
KEINJ=0.5#(0.005*WJET Jex2 3 EPINJZ0.09*KEINJH#1.5/(0.05"0.5) 
KEINA=0.5%( .005%1.0 )¥*2 3 EPINA=0.09*KEINA®*] .5/(0.1) 


GROUP 10. INTER-PHASE-TRANSFER PROCESSES AND PROPERTIES 


GROUP 11. INITIALIZATION OF VARIABLE OR POROSITY FIELDS 
FIINIT( H1 )=HAMB 
FIINIT( TMPLI=ATAMB 3 FIINIT(DENIT )J=1.0 
FIINIT(KE J=KEINA 3 FIINIT(EP JSEPINA 
FIINIT(P1)=0.0 
FIINIT(W1)=0.0 
FIINIT( V1 )=0.0 
FIINIT(U1 )=0.0 
PATCHIWINIT,INIVAL,1,5,1,8,15,30,1,1) 
INIT(WINIT,W1,0.0,0.2*WJET ) 
INIT(CWINIT,H1,0.0,CSUBPxO.2*TJET ) 
INIT(WINIT ,KE,0.0,KEINJ*0O. 2) 
INIT(WINIT,EP,0.0,EPINJ#0.2) 
wee SEE GROUP 13 FOR POROSITY FIELOS 


GROUP 12. CONVECTION ANO DIFFUSION AOJUSTMENTS 


GROUP 13. BOUNOARY CONOITIONS ANO SPECIAL SOURCES 
**#OEFINE SOURCES AND SINKS FOR JET ENGINE x* 

CONPOR( 0.0,CELL,1,1,49,5,12,19) 

PATCH( JETIN,HIGH,1,1,4,5,11,11,1,1) 

COVAL{ JETIN,P1,FIXFLU,-MJET ) 

COVALt JETIN,H1,ONLYMS,SAME ) 

PATCH( JETOUT,LOW,1,1,%,5,15,15,1,1) 

COVAL( JE TOUT,P1,FIXFLU,MJET ) 

COVAL( JE TOUT »H1,ONLYMS »HJET ) 

COVAL( JE TCUT,W1,ONLYMS ,WJET ) 

COVAL( JETOUT »KE ,ONLYMS ,KEINJ ) 

COVAL( JETOUT ,EP,ONLYMS,EPINJ } 


#*XEAST BOUNOARY TO ATMOSPHERE ¥* 
PATCH{ OUTSIOE ,EAST,10,10,1,12,1,40,1,1) 
COVAL( OUTSIOE ,P1,0.1,0.0) 

COVAL( OUTSIOE 5H1,FIXVAL,HAMB ) 


*#* TOP OF HUSH HOUSE ANO TUBE INCLUDING TUBE LIPx 
CONPOR(0.0,NORTH,1,9,12,12,1,19) 
CONPOR(0.0,NORTH,1,5,-8,-8,18,27) 

CONPOR{ 0.0,NORTH,1,5,-2,-2,35,35) 
CONPOR(0.0,NORTH>1,55-%»-4,36,36) 
CONPOR( 0.0,NORTH,1,5,-5,-5,37,37) 
CONPOR(0.0,NORTH,1,5,-°6,-6,38,38) 


- 
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CONPOR(0.0,NORTH)1,55-7)27)39)39) 


*#*TOP OF MODEL TO ATMOSPHERE ## 
PATCH( TOPOUT ,NORTH,1,10,12,12,20,40,151) 
COVAL( TOPOUT ,P1,0.1,0.0) 
COVAL( TOPOUT ,H1,ONLYMS ,HAMB ) 
PATCH( TOPP,NORTH,10,10,12,1251,19,1,1) 
COVAL(TOPP,P1,0.1,0.0) 
COVAL(TOPP ,H1,ONLYMS ,HAMB ) 


**FLOOR OF HUSH HOUSE, TUBE AND GROUNDx* 
CONPOR(0.0,SOUTH,1,10,1,-15,1,40) 


**FORE AND AFT BULKHEADS IN HUSH HOUSE ANDO TUBE + 
CONPOR(0.0,LOW,1,9,1,12,-1,-1) 
CONPOR(0.0,LOW,6,9,1,12,-20,-20) 
CONPOR(0.0,LOW,1,5,9,12,-20,-20) 
CONPOR(0.0,L0WN,1,5,1,2,-35,-35) 
CONPOR(0.0,L0OW,1,5,354%,-36)-36) 

CONPOR(0.0,LOWN,1,5,5 95,737,737) 
CONPOR(0.0,LOH,1,5,6,6,-58,-38) 
CONPOR(0.0,L0OW,1,55,7,7,7-39,-39) 
CONPOR(0.0,HIGH,1,5,8,8,-39,-39) 


*#* VERY FRONT MOOEL STRIP TO ATMOSPHERE 
PATCH( FRONT ,LOW,10,10,1,12,1,1,1,1) 
COVAL( FRONT ,P1,0.1,0.0) 
COVAL( FRONT »H1,ONLYMS , HAMB ) 


*#* THE VERY ENO OF OUR MOOEL TO ATMOSPHERE x 
PATCH( ENOENO ,HIGH,1,10,1,12,40,40,1,1) 
COVAL(ENDEND,P1,0.1,0.0) 

COVAL(ENDEND >H1,ONLYMS,HAIB ) 


*#*OUTSIOE EAST WALL IN HUSH HOUSE 
CONPOR(0.0,EAST,-9,-9,95,12,1,19) 
CONPOR(0.0,EAST,-9,-95,1,8,1,1) 
CONPOR(0.0,EAST ,-95,-951 »8,3,3) 
CONPOR(0.0,EAST,-9,-9,1,8,5,58) 

CONPOR( 0.0,EAST ,-9,-951 »8,10,11) 
CONPOR(0.0,EAST »-9,-9,1,8,13,19) 
CONPOR(0.0,EAST,-9,-95,1,8,18,19) 


**xOUTSIDE OF TUBE PLUS TUBE LIP IN HUSH HOUSE 
CONPOR(0.0,EAST,-5,-5,1,8,18,39) 


GROUP 14. DOWNSTREAM PRESSURE FOR PARAB=.TRUE. 


GROUP 15. TERMINATION OF SWEEPS 
FSWEEP=1; LSWEEP=180 


GROUP 16. TERMINATION OF ITERATIONS 


GROUP 17. UNOER-RELAXATION DEVICES 
REAL(OSLT,FACT } 
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DELT=50.0/1NZ*0.2*WJET ) 

FACT=10. 

FACT=2.5 
FACT=0.75 

FACT=3.5 

FACT=S. 

FACT=0.50 

FACT=0.25 
RELAX(P1,LINRLX,0.3) 
RELAX(U1,FALSOT,FACT*DELT ) 
RELAX(V1,FALSOT,FACT#DELT ) 
RELAX(hi1,FALSOT,FACTXDELT ) 
RELAX(H1,FALSOT,2.5*FACT*DELT ) 
RELAX(KE ,FALSOT,FACT#DELT ) 
RELAX(EP,FALSDT ,FACT*DELT ) 


GROUP 18. LIMITS ON VARIABLES OR INCREMENTS TO THEM 
VARMIN( U1] )=-300. 
VARMAX(U1)=1000. 
VARMIN( V1 )=-300. 

VARMAX(V1 )=1000. 
VARMIN(W1 )=-100. 
VARMAX(W1 )=1500. 
VARMIN( TMP 1 )=250 
VARMAXUTMP1 J=TJET 
VARMIN(H1)=1000. 
VARMAX(H1 J=CSUBP*TJET 
VARMIN( P21 )=-3.5E+05 
VARMAX(P1)= 3.5&+05 
VARMIN( RHO] )=.05 
VARMAA( RHO1 )=5. 
VARMIN(KE )=1.E-08 
VARMAX(KE )=1.3*KEJET 
VARMIN(EP )=1.E-08 
VARMIN(VIST )=0.2 
VARMAX(EP )=1.3%(KEJET##H1.5) 


GROUP 19. DATA COMMUNICATED BY SATELLITE TO GROUND 
GROUP 20. PRELIMINARY PRINT~OUT 
GROUP 21. PRINT-OUT OF VARIABLES 


GROUP 22. SPOT-VALUE PRINT-OUT 
**¥MONITOR THE OUTPUT AT DESIGNATED POINTS 
IXMON=1 3 IYMON=4 3 IZMON=16 
TSTSWP=1 3 LUPR3=6 
ABSIZ=1.0 3 ORSIZ=1.0 


GROUP 23. FIELD PRINT-OUT AND PLOT CONTROL 
NUMNCLS=10*NXPRIN 
ITAGL=3 3 NPLT=1 3 NZPRIN=l 
IPLTL=LSWEEP 
NPRINT=3 
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GROUP 2¢. DUMPS FOR RESTARTS 
SAVE =T 
RESTRTCALL) 
STOP 


cccc HHH PHOENICS VERSION 1.3, 03 SEPT 1986 


HHHH (C) COPYRIGHT 1984 

peceeee umriteitinn ES ha al baa aes MOMENTUM LTD 

CCCCC HHNHHKHHHHHH ALL RIGHTS RESER ‘ 
cee HHHHHHHHHHHHH CHAM LTD, BAKERY HOUSE, 40 HIGH ST 
CCCCCCC HHHKHHHHHHHH WIMBLEDON, LONDON, SH19 SAU 

CCCCCCC HHHHHHHHHH~ TEL: 01-947-76513 TELEX: 928517 

cccccccC HHHHH FACSIMILE: 01-879-3497 
cccc HHH THE OPTION LEVEL IS -18 


THIS CODE MAY ONLY BE USED UNDER THE TERMS AND CONDITIONS 
OF A LICENCE AGREEMENT WITH CHAM LTO. 


REPLICATION OF THIS CODE IS PROHIBITED UNLESS 


SPECIFICALLY AUTHGRISED IN WRITING BY CHAM LTD. 


He FTE IEEE TEI EIT TE FETE ITE IEE IE IEE HE FE IEE FEE EE IEEE TORINO HIE IEEE IEE 


GREX1 OF 15/07/86 HAS BEEN CALLED 
FORMATTED SATLIT DATA READ FROM OF10 FOR IRUN= 1 


KKH eonenn STORAGE INFORMATION ------ HH > 
F DIMNSN=200000 OCCUPIED=123707 ESTIMATED MINIMUM DIMNSN= 23267 
DEP VRBL= 67200 OLD VRSL= 0 30 COEFF= 23880 3D DVDPS= 14400 


INITIAL FIELOS READ FROM DFO9 

FE He IE ITE EE EE TE FETE IEF EEE IE EOI 90226 EEE TE ITE IEIE IEEE TEE EE IIE IEE IEIEIE IE 3 
GROUP 1. RUN TITLE & NUMBER 

FOR EO RE OE EIS IIE IEE IIIB 
eet eet ee eee ee ee ce eee i ee eee cd 


TEXT(NAS JACKSONVILLE HUSH CHAMBER ) 


PEI TEI TEM TIEN HE Se HE IE TE HE IERIE IE IE HE IEIETE HH IE TEE FE HE ITE ETE HE HE IEEE EEE HE EEE FE IE IEEE FEE 
BE EIT IEN Pe TE IEE HE ETE HE ETE EK HE IE TE HE IE ITE TE HEI HE IEE IEE FEI IE HE ICI FICHE HE IE HE IPE IIE NE HK FEN IE HE HE 


TRUN | = | 

HK TMK HK EH HE EI A HIE IIE HE IE TEE I FIERCE EE FO FOI III IE IE IE 5 
GROUP 2. TRANSIENCE 

STEADY = T 

MT RE IE EAR IESE HIE SIE SE te HE EEK IE RIE I RIE 2 THI IHR 
GROUP 3. X-DIRECTION GRID SPACING 


CARTES = T 

NX = 10 

XULAST = 1.000E+00 

METHOD OF PAIRS USED FOR GRID SETTING 

XFRAC ( 1) = -1.000&+00 ;XFRAC ( 2) = 3.100E-01 
XFRAC ( 3) = @.000£+00 3xXFRAC ( ©) = S5S.475E-01 
XFRAC ( 5) = 3.000E+00 ;XFRAC ( 6) = 2.000E+00 
XFRAC ( 7) = 2.000E+*00 3XFRAC ( 8) = 1.400E+00 


3 JETER SE IE IEE JE SE EEE EEE IEEE IHHEHHEE 
GROUP 4. Y-DIRECTION GRID SPACING 
NY = 12 
YVLAST = 1.000E+00 
METHOD OF PAIRS USED FOR GRID SETTING 


Wak 


YFRAC ( 1) = -3.000E+00 ;YFRAC ( 2) = 5.620E-01 
YFRAC ( 3) = 2.000E+00 sYFRAC ( 4%) = 3.140E-01 
YFRAC ( 5) = 3.000E+00 s3YFRAC ( 6) = 5.620E-01 
YFRAC ( 7) = G.QOOE+00 3YFRAC ( 8) = 7.500E-01 


OOO IR a ICI IC IIIB On RII FEIHIOR BIBI IRIS CRORE He 
GROUP 5. Z-DIRECTION GRID SPACING 


PARAB = F 
NZ = 40 
ZHLAST = 1.000£+00 


METHOD OF PAIRS USED FOR GRID SETTING 


ZFRAC ( 1) = ©1.000E+00 3Z2FRAC ( 2) © 1.000E+00 
ZFRAC ( 3) = 1.Q00E+00 3ZFRAC ( 4) © 1.700E+00 
ZFRAC ( 5) = 1.000E+00 3ZFRAC ( 6) = 2.000E+00 
ZFRAC ( 7) = 1.000E+00 3ZFRAC ( 8) = 1.700E+00 
ZFRAC ( 9) = 2.000E+00 3;Z2FRAC {( 10) = 3.000E£+00 
ZFRAC ( 11) = 1.Q000E+00 3ZFRAC ( 12) = 1.500E+00 
ZFRAC ( 13) = J1.000E+00 3;ZFRAC ( 14) = 1.000E+00 
ZFRAC ( 15) = 12.000E+00 3ZFRAC ( 16) = 1.700E+00 
ZFRAC ( 17) = 2.000E+00 3ZFRAC ( 18) = 1.000E+CO 
ZFRAC ( 19) = J1.000E+00 ;ZFRAC ( 20) = 1.700E+00 
ZFRAC ( 21) = 2.000E+00 3ZFRAC ( 22) = 1.000E+00 
ZFRAC ( 23) = 3.Q000E+00 3ZFRAC ( 24) = 5.667E-01 
ZFRAC ( 25) = 2.000E+00 3ZFRAC ( 26) = 3.500E-01 
ZFRAC ( 27) = 7.Q000E+00 sZFRAC ( 28) = 3.000E+00 
ZFRAC ( 29) = 1.000E+00 3ZFRAC ( 30) = 2.000E+00 
ZFRAC ( 31) = G.OQOOQO0E+00 3ZFRAC ( 32) = 7.500E-01 
ZFRAC ( 33) = %.OQO0E+00 3ZFRAC ( 3%) = 5§.750E-01 
ZFRAC ( 35) = 4.000E+00 3ZFRAC ( 36) = 5§.Q00E-01 
ZFRAC ( 37) = 1.000E+00 s3ZFRAC ( 38) = 3.000E+00 


HE HHH HH WH HH He HI HIE Th HE HE Te eH eH HE HEI ITE EI IE EI HE EIEN IE IEEE EEE 
GROUP 6. BFC AND GEOMETRY CHANGES 

MH HH HIM ITM He PEE EHH HIM HIE HM HH HH HIE 5k IEE IOI HH AIH IE FE HIE HE HE HEHEHE HEI HE 
GROUP 7. VARIABLES: STOREO,SOLVEO,NAMED 

ONEPHS = 7 

NAME( 1) =P1 

NANE( 5) =V1 


sNAME( 3) =U1 
sNANE( 7) =H1 
NAME( 12) =KE psNAME(13) =EP 
NAME(1¢) SHI] sNAME(19) =HPOR 
NAME( 20) =NPCR ;NANE( 21) =EPOR 
NANE( 22) SVPOR 3NAME(23) =VIST 


NAME(24) =RHOL 


SOLUTN(P1 
SOLUTN(U1 
SOLUTN(V1 
SOLUTN(W1 
SOLUTN(KE 
SOLUTN(CEP 
SOLUTN(H1 


»¥ YY N,N,N) 
2V¥sY,N N,N,N) 
»>Y»VY>N>»N»N,N) 
»¥ »Y»N»N N,N) 
»Y»Y,N,N»N,N) 
»Y¥»Y»N N,N,N) 
»Y¥»¥ »N N,N,N) 


SOLUTN(HPOR,Y»N»N»N»N JN) 
SOLUTN(NPOR,Y,N,N»N,N>N) 
SOLUTN(EPOR,Y>»N),N, N,N,N) 
SOLUTN( VPOR>YsN,»N»N NON) 
SOLUTNI VIST,>Y,N»N N N,N) 
SOLUTN( RHO1>Y>N»N»N NIN) 
SOLUTN( TMP1,Y N,N NN) 


DENI 
VIST 
EPOR 
TEMP 


2% 
23 
21 35HPOR 
25 


sNAME( 25) =TMP1 


19 5NPOR = 20 3VPOR = ce 


FRR REHASH ARR 
GROUP 8. TERMS & DEVICES 


.OOOE-O01 ;ZDIFAC = 1.000E+00 


F 


HE HE He He IIE HE HE HE HE IE IE IE IE IE FE FE IE IE IEIIHIE FHIHIHIIEIHIIOIHIE NE 2 IIE IEE IEE IE FE HE FE IE HE HE HEE HIE 


1 = 


1B CO 


TERMS (Pl ,Y¥r¥,V»N»Y>N) 
TERMS (Ul ,Y,Y¥>¥sN>Y,N) 
TERMS (V1 5Y,YsY»N>Y2N) 
TERMS (HL ,YsYsYoN,rY>N) 
TERMS (KE aN r>VrVoNrYN) 
TERMS (EP 5 N»rY,Y>N,YN) 
TERMS (Hl »Yr¥>Y»N>Y>N) 
DIFCUT = 5 
GALA = F j;AOQDOIF = 
NEWRH1L = T 

NEWENT = T 

GROUP 9. PROPERTIES 
RHO1 =-1.016E+04 3;TMP 
RHOI1A = 0.000£+00 ;RHO1B = 
PRESSO™ = ]-027S5E+C5 
TMP1A = 1.000E-20 3TMP 
TEMPO = 0.000E+00 

ELIA = 0.000E+00 ,EL1B 


-1.013&+0% 3ELI =-1.015E+04 
3.496E-03 3;RHOIC = 7.143€-01 
9.960E-0% ;TMP1C = 0.000&+00 
0.000E+00 ;ELIC = 0.000E+00 
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ENUL = OOOE-05 ;ENUT =-1.014E+04 

ENUTA = OOO0E+00 ;ENUTB = 0.000E+00 ;ENUTC 
DRH1DP =-1.016E +04 

PRNOTL(U1 = 1.000&+00 3;PRNDTL(V1 } = 1.000£+00 
PRNOTLIWI = 1.000E+00 ;PRNDOTL(KE ) = 1.000&+00 
PRNDOTLCEP = 1.000E+00 ;PRNOTL(H1 ) = 1.000£&+00 
PRT (U1 = 1.000&+00 3PRT (Vl )} = 1.000&+00 
PRT (W1 = 1.00CE+00 3PRT (KE J) = 1.000E+00 
PRT (EP = 1.314E&+00 3PRT (Hl ) = 1.000E+00 


= 0.000E+00 


ERE IR RE poe ETE EEOC IIE DE SEBS EI IRIE aE 
GROUP 1O.INTER-PHASE TRANSFER PROCESSES 


FE FE IE IEE IE IEE HEH HE HIE IIE EEE SIE FE HEE FETE FE ETE IEE IEE IE HE HE IE FE IIE FE IONE IEE HII 


GROUP LTL.INITIALIZE VAR/POROSITY FIELDS 


FIINIT(P1L 
FIINIT(V1 
FIINIT(KE 
FIINIT(H1 


FIINIT(NPOR ) 
FIINIT(VPOR ) 
FIINIT( RHO] ) 


PATCH(SOUTH1 
INIT( SOUTH 
INIT(SOUTH1 


PATCH( LOW] 
INIT( LOW 
INITCLONY 


PATCHICELL1 


PNITECELEL 
INIT(CELL1 
PINE RCCELLT 
INIT(CELLI 
INIT(CELLI 
INIT(CELLi 
INIT(CELLI 
INIT(CELLI 
INIT(CELLI 
THITCCELLY 


= -1.022£ +04 
-1.022E +04 
-1.022E +04 
-1.022E +04 
-1.022E +04 
-1.022E +04 
-1.022E+04% 


te | | | 


»INIVAL » 1; 
» 0.000E+00, 0.000E+00 ) 
»NPOR, 0.000E+00, 0.000E+00) 


»V1 


»LNIVAL > l; 
» 0.000E&+00, 
»HPOR, 0.000E+00, 


»W1 


»INIVAL> l> 


»Pl 
»U1 
»V1 
»R1 
oKE 
»EP 
»H1 


? 


wv w we we Ww 


? 


»HPOR, 
»NPOR >» 
»EPOR, 


0.000E+00, 
0.000E+00, 
0.000E+00, 
0.000E+00, 
0.000£+00, 
0.000E+00, 
0.000E +00, 
0.000E+00, 
0.000E+00, 
0.000E+00, 


sFIINIT(UlL § ) 
SFIINIT(WH1 = ) 
SFIINIT(EP § ) 
sFIINITCHPOR ) 
sFIINITCEPOR ) 
SFIINITC(VIST ) 
sFIINITC(TMP1 ) 


1; 3, 3> 


i; G, 5; 


0.000E+00 ) 
0.000E +00) 


1, %» 5, 


0.000E+00 ) 
0.000E +00) 
0.000E+00 ) 
0.000E+00 ) 
0.000E+00 ) 
0.000£ +00) 
0.000E+00 } 
0.000£+00 ) 
0.000E +00 ) 
0.000E +00 ) 
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-1.022E +04 
-1.022E +04 
-1.022E +04 
-1.022E+04¢ 
-1.022E +04 
-1.022E£ +04 
-1.022E+04 


12, 14, 


ll, 


1l; 


12, 14; 


1, 


1; 


i; 


1) 


1) 


7) 


PATCH(EAST21 ,INIVAL, 9» 9» 1> 8, 3> 3, 1; 1) 
INIT(EAST2Z1 ,U1l , 0.000E+00, 0.000E+00) 
INIT(CEAST21 ,EPOR, 0.000E+00, 0.000E+00) 


PATCH(EAST22 ,INIVAL, 9, 9, l> 8, 5S» 8, 1> 1) 
INIT(EAST22 ,U1l , 0.000E+00, 0.000E+00) 
INIT(EAST22 ,EPOR, 0.000E+00, 0.000E+00 ) 


PATCH(EAST23 ,INIVAL; 9» 9» l; 8, 10, ll, 1, 1) 
INIT(EAST23 ,U1l , 0.000E+00, 0.000E+00) 
INIT(EAST23 ,EPOR, 0.000E+00, 0.000E+00) 


PATCH(EAST2¢@ ,INIVAL, 9» 9» l> 8, 13, 14, 1; 1) 
INIT(EAST2%¢ ,U1l , 0.000E+00, 0.000E+00) 
INIT(EAST2% ,EPOR, 0.000E+00, 0.000E +00) 


PATCH(EAST2S5 ,INIVAL, 9, 9, l> 8 18, 19, 1l> 1) 
INIT(EAST2Z5 ,U1l , O.000E+00, 0.000E+00 ) 
INIT(EAST25 ,EPOR, 0.000E+00, 0.000E+00) 


PATCH(EAST26 ,INIVAL, 5, §, 1, 8, 18, 39, 1, 1) 
INIT(EAST26 ,U1l , 0.000E+00, 0.000E+00) 

INITCEASTZ2Z6 ,EPOR, 0.000E+00, 0.000E+00 ) 

INIADD = F 

LUFI = 9 : 

NAMFI =DFO9 

FEE II EI OEE ICRI RI CSE LTE FIFI IRIHIEIE IH REE 
GROUP 13. BOUNDARY & SPECIAL SOURCES 


PATCH(KESOURCE ,PHASEM, 1, 10, 1, 12> 1, 40, i; 1) 
COVALIKESOURCE ,KE ,-1.015E+04,-1.015E +04) 
COVAL(KESOURCE,EP »-1.015E+04,-1.015E +04 ) 


PATCH( JETIN »HIGH 1; 1» 4) §, ll, 11, i> 1) 
COVAL(JETIN »PL1 » 1.000E-10,-1.977E+02) 
COVAL( JETIN vsH1l , 0.000E+00,-1.026E +04) 


ee w 


PATCH(JETOUT ,LOW l; 1; 4, 5» 15, 15, 1, 1) 
COVAL(JETOUT ,P1l » 1.Q00E-10, 1.977E+02) 


w 


Us 


COVAL(HWALL12 
COVAL(HWALL12 
COVAL(HWALL12 
COVAL(HWALL12 


PATCH(LRALL12 
COVAL(LHALL12 
COVAL(LWALL12 
COVAL(LWALL12 
COVAL(LWALL12 


PATCH( HWALL13 
COVAL(HWALL13 
COVAL( HWALL13 
COVAL(HRALL13 
COVAL( HWALL13 


PATCH(LWALL13 
COVAL(LWALL13 
COVAL(LWALL13 
COVAL(LWALL1I3 
COVAL( LWALL13 


PATCH( HWALL14 
COVAL(HNALL14 
COVAL(HNALL14 
COVAL(HHALL14 
COVALIHMALL14 


PATCH(LWALL14% 
COVAL(LWALL14 
COVALILWALL1¢ 
COVAL(LWALL1¢ 
COVAL(LWALL14 


PATCH(HWALL15 
COVAL(HWNALL15 
COVAL(HWALLI5 
COVALCUHWALLI1S 
COVAL(UENALLI5 


2»U1l »-1.013E+04¢, 0.000E+00 ) 
»V1l_s- 9 -1.013E +04, 0.000E +00 ) 
»9KE »-1.012E+04¢,-1.013E +04) 
»9EP = 9-1.013E404,-1.013E +04) 


»LWALL ;» 


5» 9, 12, 


»Ul »°1.013E+04¢, 0.000E+00 ) 
»Vl 971.013E+04, 0.000E+00 ) 
KE 971.013E +04 »-1,.013E +04) 
»9EP) 421.013E+04,-1.015E +04) 


gHWALL 5 


5> l, 2> 


»U1l = »7-1.013E+04¢, 0.000E +00 ) 
»V1_s- 45-1 .013E 404%, 0.000E +00) 
»9KE »971.013E+04,-1.013E+04%) 
»EP 971.013E +04 ,-1.013E +04) 


»LWALL > 


5, 1> 2> 


e>UlL = »7-1.013E+04, 0.000E+00) 
»Vl_s-9 -1.013E +04, 0.000E +00 ) 
»9KE 52°1.013E4+04,-1.013E+04) 
»9EP 5 -1.013E#04,-°1.013E +04 ) 


»sHWALL > 


5» 3, 4%» 


yUl  »9-1.013E+04¢, 0,.000E+00 ) 
»Vl_»21.013E+04%, 0.000E +00) 
x9KE 3 »-1.013E +04 ,-1.013E +0%) 
y9EP y= 1 .013E4+04,-1.013E +04) 


gLWALL > 


5, 3» G, 


»Ul »-1.013E+04¢, 0.000E+00) 
»Vl »°1.013E+0¢, 0.000E +00 ) 
»9KE 3 47-1.013E404¢,-1.013E +04) 
v9EP = 921.013E4+04,-1.013E +04 ) 


sHWALL 5 


5S» 5, 5, 


21 »-1.013E+*0G¢, 0.000E+00 ) 


»Vl 4» -1.013E+04¢, 0.000E+00 ) 


x9KE 3 9271.013E404,-1.013E +0% ) 
y9EP 1 .013E404,-1.015E +04 ) 
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20, 


34% 


35> 


35> 


36> 


36, 


20; 


34> 


35> 


35> 


36, 


1; 


1; 


1, 


1; 


1; 


1; 


1) 


1) 


1) 


1) 


13 


13 


PATCH(LWALL15 
COVAL(LWALLI5 
COVAL(LWALLIS 
COVAL(LWALLI5 
COVAL(LWALL15 


PATCH( HWALL16 
COVAL(HHALL16 
COVAL(HWALL16 
COVAL(HHALL16 
COVAL(HWALL16 


PATCH(LWALL16 
COVAL(LWALL16 
COVAL(LWALL16 
COVAL(LWALL16 
COVAL(LWALL16 


PATCH(HWALL17 
COVAL(HWALL17 
COVAL(HWALL17 
COVAL(HWALL17 
COVAL(HWALL17 


PATCH( LWALL17 
COVALILWALL17 
COVAL(LWALL17 
COVAL(LWALL17 
COVAL(LWALL17 


PATCH(HWALL18 
COVALCHWALL13 
COVAL(HWALLI8 
COVAL(HWALL18 
COVAL(HRALL18 


PATCH(LWALLI8 
COVAL(LWALL18 
COVAL(LHALL18 


yLWALL ,» l; 5» IF 5; 


»Ul 
»V1 


»KE 


»EP 


»71.013E +04, 0.000E+00 ) 
»71.013E+04¢, 0.000E+00 ) 
971.013E+04,-1.013E +04) 
971.013E+04,-1.013E +04) 


»HWALL > 1» 5» 6» 6> 


Ul 
»V1 


»KE 


»EP 


971.013E+0G¢, 0.000E+00) 
971.013E+04, 0.000E+00) 
971.013E+04,-1.013E£+04) 
971.013E+04,-1.013E +04) 


»LAWALL 5) l> 5, 6, 6» 


2»Ul 
»V1 
»KE 
»EP 


y71.013E+04, 0.000E+00 ) 
971.013E+04¢, 0.000E+00) 
971.013E +04 ,-1.013E+04 ) 
»7-1.013E+04,-1.012E +04) 


»aHHALL > l> 5; 7» 7> 


»Ul 
»V1 
oKE 
»EP 


»71.013E+0¢, 0.000E+00) 
»-1.013E+04¢, 0.000E+00 ) 
971.013E+04,-1.013E+04) 
971.013E*04,-1.013E +04) 


»LWHALL , 1; 5» 7> 7> 


2U] 
V1 
»KE 
»EP 


»7l.013E*04, 0.000E +00) 
»71.013E+0%, 0.000E+00 ) 
»71.013E+04,-1.013E +04) 
971.013E+04,-1 .0135E +04) 


»sHWALL , i 5> 8, 8, 


»U1 
»V1 


»KE 


»EP 


»7-1.013E+04, 0.000E+00 ) 
971.013E+04%, 0.000E+00) 
971.013E +04, -1.013E+0%) 
971.013E+04¢,-1.013E +04) 


»LWALL , l> 5; 8, 8, 


Ul 


»V1 


»71.013&E+0¢, 0.000E +00) 
971 .013E+0%, 0.000E+00 ) 


ue 


37, 


37, 


38, 


38, 


39> 


39; 


G0; 


37, 


37> 


38, 


38, 


39, 


39; 


40, 


1; 


1» 


1» 


1» 


i> 


1; 


1; 


1) 


1) 


1) 


1) 


1) 


1) 


1) 


COVAL(LWALLI8 ,KE = »-1.013£+04,-1.013E +04} 
COVAL(LWALLI8 ,EP ,-1.013E+04,-1.013E +04) 


PATCH( FRONT » LOW » 10, 10; 1, 12, 1; 1; 
COVAL( FRONT »Pl » 1.000E-01, 0.000E+00 ) 
COVAL( FRONT »HL », 0.000E+00, 2.962E+05) 


PATCH(ENDEND ,HIGH =, 1, 10, 1, 12, 40, 40, 
COVAL(ENDEND ,Pl , 1.000E-01, 0.000E+00) 
COVAL(ENDEND ,H1l », O.000E*00, 2.962E+05) 


PATCH(EWALLI19 ,ERALL , 9» 9, 9, 12; 1, 19; 
COVAL(EWALL19 ,V1l ,-1.013E+04, 0.000E+00) 
COVAL(EWALLI9 ,W1l = ,-1.013E+04, 0.000E+00 ) 
COVAL(EWALL19 ,KE »-1.013E+04,-1.013E +04) 
COVAL(EWALL19 ,EP 5-1.013E +04,-1.013E+04) 


PATCH(WHALL19 ,WHALL ,», 10, 10, 9, 12, 1, 19; 
COVALIWHALL19 »,V1l ,-1.013E+04, 0.000E+00) 
COVALCWNALL19O ,W1l = »5-1.013E+04%, 0.000E+0C ) 
COVAL(WNALLIO ,KE  ,-1.013E+04%,-1.013E +04) 
COVALEWHALLIO ,EP  »-1.013E404,-1.013E +04) 


PATCH(EWALL20 ,EWALL , 9, 9» 1; 8, 1, 1; 
COVAL(EWALL20 ,V1l ,-1.013E£+04, 0.000E+00) 
COVAL(EWALL20 »,W1l 4,-1.013E+04%, 0.000E+00) 
COVALCEWALL20 ,KE »-1.013E£+04,-1.013E +04) 
COVAL(ENALL20 ,EP ,-1.013E£+04,-1.013E +04) 


PATCH(WHALL20 ,WHWALL », 10, 10, 1; 8, 1; 1; 
COVAL(WHALL2O ,V1l ,-1.013£+0%, 0.000E+00} 
COVAL(WHALL20 »,W1l ,-1.015E+04, 0.000E+00) 
COVALIWHALL20 ,»KE »,-1.013E+09,-1.013£ +04) 
COVAL(WHALL20 ,EP 3 ,-1.013E+04,-1.013E +04) 


PATCH(EWALL21 ,EWALL , 9, 9, 1; 8, 3» 3, 
COVAL(EWALL21 ,V1l ,-1.013£+04, 0.000E+00) 
COVAL(EWALL21 ,W1l 5,-1.013E+0%, 0.000E+00) 
COVAL(EHALL21 ,KE »,-1.013E +04 ,-1.013E +04) 
COVAL(ENALL21 ,EP 4-1.013&+04,-1.015E+04) 
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PATCH( WHALL21 
COVAL(WAALL21 
COVAL(WHALL21 
COVAL(WHALL21 
COVAL(WWALL21 


PATCH(EWALL22 
COVAL(EWALL22 
COVAL(EWALL22 
COVALI(EWALL22 
COVAL(EWALL22 


PATCH(WHALL22 
COVAL(WHALL22 
COVALIWWALL22 
COVAL(WRALL22 
COVAL( WHALL22 


PATCH( EWALL23 
COVAL(EWALL23 
COVAL(EWALL23 
COVAL( EWALL23 
COVAL(ERALL23 


PATCHIWRALL23 
COVAL( WWALL23 
COVAL(WHALL23 
COVAL( WHALL25 
COVAL(WAYALL23 


PATCHIEWALL24 
COVALIEWALL24¢ 
COVAL(EWALL24 
COVALIEWALL24¢ 
COVAL( EWALL24 


PATCHIUINWALL2¢ 
COVALIWNALL24 
COVALIFWALL24¢ 


»WWALL ,» 10, 10, GS 8; 
»V1l4»-1.013E+0¢4¢, 0.000E+00) 
»N1l = »-1.013E+0%, 0.000E+00) 
»9KE »-1.013E+04,-1.013E +04) 
»EP 971.013E +04,-1.013E +04) 


»EWALL > 9; 9, l> 8, 


»V1l »-1.013E+04, 0.000E+00) 
yAl ,°1.013E+04¢, 0.000E+00) 
»9KE 57-1.013E4+04,-1.013E +04 ) 
»9EP = 9 $1 .013E404,-1.013E +04) 
»AWALL > 10; 10, diy 8; 
xl ,»-1.013E+0%, 0.000E+00) 
»yAl = ,-1.013E+04¢, 0.000E+00) 
»9KE 5-1.013E*04,-1.013E404) 
»9EP =» 2-1.013E+04,-1.013E +04) 


»EWALL >» 9, 9» 1; 8, 


»V1l_s 9 -1.013E£+04¢, 0.000E+00) 
»Al = »-1.013E+04, 0.000E+00) 
o9KE »71.013E404,-1.013E +04 ) 
y9EP =» -1 .013E +04,-1.013E +04) 
»AWALL ,» 10, 10; 1; 8, 
»V1l 4-1.013E+04, 0.000E+00) 
yAl = 4°1.0135+04, 0.000E+00) 
»9KE »°-1.013E4+04,-1.013E +04) 
»EP 4-1.013E404,-1.013E 46%) 


2»EWALL ,» 9, 9> 1; 8, 


»V1l = »-1.013E+04¢, 0.000E+00) 
»Al = »-1.013E+04, 0.000E+00 ) 
»9KE 45-1.013E404,-1,.013E +04) 
»9EP = 9 -1.013E+04¢,-1.013E +04 ) 
»NNALL » 10, 10; l> 8, 
»V1l4-1.013E+0G, 0.000E+00 ) 
yH1l 5°1.013E+04, 0.000E+00 ) 


Jez 


3, 


5, 


5; 


10, 


10; 


13, 


13; 


3» 


8, 


8, 


ll, 


ll, 


14; 


14, 


1; 


1» 


1» 


1; 


1; 


i; 


1; 


1) 


1) 


1) 


i 


1) 


1) 


1) 


COVALIWHALL2% »,KE ,-1.013E+04,-1.013E+04) 

COVAL(WHALL24 ,EP ,-1.013E+04,-1.013E +04) 

PATCH(EWALL25S ,EWALL , 95 9, 1; 8, 18, 19, 1, 
COVAL(EWALL25 »V1l  ,-1.013E+04, 0.000E+00) 

COVAL(EMALL2S »W1l 5-1.013E+04, 0.000E+00) 

COVAL(EWALL2S ,KE ,-1.013E404,-1.013E +04) 

COVAL(EWALL25S ,EP 971.013E+04,-1.013E +04) 

PATCH(WHALL2S ,WNALL », 10, 10, 1; 8, 18, 19, 1, 
COVALUWWALL2ZS »,V1l ,-1.013E+04¢, 0.000E +00) 

COVAL(UWNALL2S ,»W1l = ,-1.013E+0%, 0.000E +00) 

COVAL(WWALL2S ,KE ,-1.013E+04,-1.013E +04) 

COVAL(WHALL2S ,EP ,-2.013E+04,-1.013E +04) 

PATCH(EWALL26 ,EWALL , 5, 5, l; 8, 18, 39, 1; 
COVAL(EWALL26 ,Vl ,-1.013E+04, 0.000E+00) 

COVAL(EWALL26 ,W1 ,-1.013E+0%, 0.000E +00) 

COVAL(EWALL26 ,»KE »,-1.013E+0¢,-1.013E +04) 

COVAL(EWALL26 ,EP ,-1.013E+04,-1.013E +04) 

PATCH(WWALL26 ,WWALL , 6, 6> 1, 8» 18), 39; i> 
COVAL(IWWALL26 ,V1l  ,-1.013E+0¢, 0.000E+00) 

COVALIWWALL26 ,»H1l ,-1.013E+04, 0.000E+00) 

COVAL(WWALL26 KE 4-1.013E+04¢,-1.013E+04) 

COVALUWWALL26 ,EP ,-1.013E+04,-1.013E +04) 

XCYCLE = F 

KELIN = 8) 


FEI HEI HE HE IE FEM FE HE ETE HE HE FETE HE IEE ITE HE IEICE 1b IEE IEICE HE FEE IE FI III FOIE 


GROUP 14. DOWNSTREAM PRESSURE FOR PARAB 
96 3 HH EET IE EE FEE EE EE EEE JHE EEE III FEI III FOIE IEE 


GROUP 15. TERMINATE SWEEPS 

LSWEEP = 180 ;ISWCl = 1 

LITHYD = 1 5LITFLX = 1 >LITC = 1 3ITHC1l = 
RESREF(P1l ) = 1.Q000E-06 3RESREF(UL ) = 1.000E-06 
RESREF(V1l ) = 1.000E-06 ;RESREF(W1 ) = 2.000E-06 
RESREF(KE ) = 1.Q000E-06 ;RESREF(EP ) = 1.000E-06 
RESREF(H1l ) = 1.Q000E-06 


FE HH EH He HHH HH HII HE IE IEE HE EIEIO IIE EI IE KIRIN HH HII HK HK ION 


GROUP 16. TERMINATE ITERATIONS 
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LITER (Pl ) = 2O ;LITER (UL ) = 20 

LITER (Vl ) = 20 5LITER (Wl ) = 20 

LITER (KE ) = 20 3LITER (EP ) = 20 

LITER (H1 ) = 20 

ENDIT (Pl ) = 1.000E-06 ;ENOIT (Ul J) = 1.000E-06 
ENDIT (Vl ) = J.000E-06 ;ENOIT (Wl ) = 1.000E-06 
ENODIT (KE ) = 1.Q000E-06 ;ENCIT (EP ) = 1.000E-06 
ENDIT (Hl ) = 1.000E-06 


FE RR IERIE FIFI IEE EIEIO RRP 
GROUP 17. RELAXATION 


RELAX(P1 ,LINRLX, 3.Q000E-01) 
RELAX(U1l ,FALSDT, %.321E-03) 
RELAX(V1l ,FALSDT, 4.321E-03) 
RELAX(WL ,FALSOT, 4.321E-03) 
RELAX(KE ,FALSDT, 4.321E-03) 
RELAX(EP ,FALSOT, %.321E-03) 
RELAX(H1 ,FALSOT, 1.080E-02) 


RELAX(VIST,FALSOT, 
RELAX(RHO1L,FALSOT, 1.000E+00 ) 
RELAX(TMP1,FALSOT, 1.000E+00) 
FEI Se Fe OR SE HO RO IOI IO IGIS SII IIB OIE HEHEHE GBH tr 


1.Q000E +00) 


GROUP 18. LIMITS 
VARMAX(P1l ) = 3.500E+05 ;VARMIN(P1 ) =-3.500E+05 
VARMAX(Ul ) = 1.000E+03 3VARMIN(UL J) 2-3.000E+02 
VARMAX(V1l ) = 1.Q000E+¢03 3sVARMIN(V1 ) =-3.000E +02 
VARMAX(W1l ) = 1.500&+03 3$VARMIN(H1 3 =-1.000E+02 
VARMAX(KE ) = 1.Q000E4¢10 ;VARMIN(KE ) = 1.000E-10 
VARMAX(EP ) = 1.Q0Q00E+10 3;VARMIN(EP ) = 1.Q0Q00E-10 
VARMAX(H1l ) = 1.952E+*06 3VARMIN(H1 ) = 1.000E+03 
VARMAX(HPOR) = 1.Q000E+10 ;VARMIN(HPOR) =-1.000E+10 
VARMAX(NPOR) = 1.000E+10 3VARMIN(NPOR) =-1.000E+10 
VARMAX(EPOR}) = 1.Q000E+10 3;VARMIN(EPOR) =-1.000E+10 
VARMAX(VPOR) = 1.000E+410 s;VARMIN(VPOR) =-1.000E+10 
VARMAX(VIST} = 1.000E4¢10 sVARMIN(VIST) = 1.Q00E-20 
VARMAX(RHO1) = 5.Q00E+00 sVARMIN(RHO1L) = 5.000E-02 
VARMAX(TMP1) = 1.94¢E+03 sVARMIN(TMP1) = 2.500E+02 


SRR REO REE III FRET IL SITE SERIE Ra J 
GROUP 19. EARTH CALLS TO GROUND STATION 

USEGRD = T j;USEGRX = T 

NAMGROD =NONE 


Gi 


GEIK = =—tC«<~ST 

FE 3 IE IE IIE EE ESE 3 IOI FE FBI RRIF JRE ISHN 99 IEE EE 
GROUP 20. PRELIMINARY PRINTOUT 

ECHO = fT 

PERI IE 99 IIE IIE ICE 9638 FIFE 36 EE 3b FEES IIE IEE 2 IEF IE FEE HE IEE FETE EDL FE FEE IEEE IE 
GROUP 21. PRINT-OUT OF VARIABLES 

INIFLO = FF j3SUBWGR= fF 

OUTPUT(P1 5YsNyNoYoY>¥3 

OUTPUT(UL ,YsNoNoYsY,¥) 

OUTPUT(V1l »YsNoNoYoYoY) 

OUTPUTIWN1 sYsNoNoYoYoY) 


OUTPUTIKE »Y,NoN,Y,Y>Y) 
OUTPUT(EP ,Y,NN,Y,Y,Y) 
OUTPUT(H1l ,Y,NN,Y,Y,Y) 
OUTPUT(HPOR,N,N N,N,N IN) 
OUTPUT (NPOR>,N»N»N>N,NDN) 
OUTPUT(EPOR,N,N»N»N NN) 
OUTPUT(VPOR,N»N»N,N»N,N) 
OUTPUT(VIST,Y,N»N»N»NIN) 
OUTPUT (RHO1,Y »N»N»N NN) 
OUTPUT(TMP1,Y>,N,N,N NIN) 


FORO ROR RS i IGE IS III UE IIB GOI IRI REE 
GROUP 22. MONITOR PRINT -OUT 


IXMON = 2 3;IYMON = G 3IZMON = 16 
NPRMON = 10000 3NPRMNT = 10000 3TSTSWP = 1 
HIGHLO = F 


U Saiekotet dk katsnedel- iodedat din -tkensdoded tok -acksdehehakedehakeh hekehedehideketehtotehstateteds aed takakehoded 
GROUP 23.FIELD PRINT-OUT & PLOT CONTROL 


NPRINT = 10000 ;NUMCLS = 10 

NXPRIN = 1 ;IXPRF = 1 sIXPRL = 10 
NYPRIN = 21 3IYPRF = 1 ;IYPRL = 12 
NZPRIN = 1 ;IZPRF = 1 3IZPRL = 40 
XZPR = F 5YZPR = F 

IPLTF = 1 ;IPLTL = 180 3sNPLT = 1 
ITABL = 3 3IPROF = 1 

ABSIZ = 1,.000E+00 sORSIZ = 1.000&+00 

NTZPRF = 1 sNCOLPF = 50 

ICHR = 2 sNCOLCO = 45 s;NROWCO = 20 


NO PATCHES YET USEQ FOR THIS GROUP 
FE IE II He I I IEE IE IEE HE EEE IEE GE IE IE ETE FE ETE TE ESE IEE TE IE TES 9 IEE 3638 3 IEE 
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GROUP 24. DUMPS FOR RESTARTS 
SAVE = T  jAUTOPS = F 
NSAVE =DFO9 
LUSAVE = “4 


sNOWIPE = F 


we GRID-GEOMETRY INFORMATION si 
RX-COOROINATES OF THE CELL CENTRES 


Te 550E-0] 985.8376 =01 1.2131£+00 1.679& +00 2.226E +00 
3.500£& +00 5.500E+00 7.500E+00 9.200E +00 1.060E+01 

Y-COOROINATES OF THE CELL CENTRES 
2.810E-01 8.430E-01 1.%05E +00 1.843E +00 2.157E +00 
2.595E+00 3.157E +00 3.719E&+00 4.375&+00 5.125E& +00 
5.875£&+00 6.625£&+00 

Z-COORDINATES OF THE CELL CENTRES 
5.000E-01 1.850E +00 3.700E&+00 § .550E +00 7.900E+00 
1.090E +01 1.315E +01 1.440E +01 1.575E€ +01 1.710E+01 
1.810E+01 1.945&+01 2.080E +01 2.180E+01 2.2586+01 
2.315E +01 2.372E+01 2.917E +01 2.452E +01 2.620E+01 
2.920E+01 3.220E+01 3.520&+01 3.820£+01 4.120E+01 
4.420E+01 4.670E+01 4.897E+01 4.882E+01 4.957E +01 
§.032E+01 §.099&+01 5§.156E+01 §.214E+01 §.271E&+01 
§.325€+01 §.375€+01 5§.4255+01 §.475E+01 5.650E +01 

~-- INTEGRATION OF EQUATIONS BEGINS --- 

TIME STEP = 1 SHEEP = 2 

TOTAL RESIOUAL/( 1.000£-06) FOR Pl IS 8.007E+05 

TOTAL RESIOUAL/( 1.000E-06) FOR Ul IS 9.092E+06 

TOTAL RESIDUAL/( 1.000E-06) FOR V1 IS 1.282&+07 

TOTAL RESIOUVAL/( 1.000E-06) FOR H1 IS 2.376&+07 

TOTAL RESIOUAL/( 1.000E-06) FOR KE IS 1.210&+09 

TOTAL RESIOUAL/( 1.000E-06) FOR EP IS 1.490E411 

TOTAL RESIOUAL/( 1.000E-06) FOR H1 IS 5.642E+10 

TIME STEP = 2 SHEEP = 3 

TOTAL RESIOUAL/( 1.000E-06) FOR Pl IS 3.733E+405 

TOTAL RESIOUAL/( 1.000E-06) FOR Ul IS 1.050E+07 

TOTAL RESIDUAL/( 1.000E-06) FCR V1 IS 1.4%74E+07 

TOTAL RESIDUAL/( 1.000E-06) FOR Wl IS 2.037E+07 

TOTAL RESIDUAL/( 1.0COE-06) FOR KE IS 1.424E+09 


BS 


TOTAL RESIDUAL/( 
TOTAL RESIOUALZ( 
TIME STEP = 1 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIDUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TIME STEP = a 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUALZ( 
TOTAL RESIOUAL/( 
TOTAL RESIOUALZ( 
TOTAL RESIOUAL/( 
TIME STEP = 1 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUALZ( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TIME STEP = 1 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUALZ( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 
TIME STEP = i 
TOTAL RESIOQUAL/( 
TOTAL RESICUALZ( 
TOTAL RESIOQUALZ( 
TOTAL RESIOUAL/( 
TOTAL RESIOUAL/( 


a ood oot oot oo oo 


ais 
1% 


1. 
1. 
1. 
1. 
1. 
1. 
1. 


1. 
.000E-06 ) 
.O00E-06 ) 
-000E-06) 
.O00E-06 ) 
-O00E-06 3 
- 0OO0E-06 ) 


eo oe Oe ee 


ee ee oe oe 


ee eo 


OO0E-06 ) 
OO0E-06 ) 
SWEEP = 
OO0E -06 ) 
O0O0E-06 ) 
O00E-06) 
OO0E-06 } 
OO0E-06 ) 
OO00E-06}) 
OO0E-06 ) 
SWEEP = 
O00E-06 ) 


SWEEP = 


.O00E-06 ) 
-O00E-06 ) 
-O00E-06 ) 
-O00E-06 ) 
.000E-06 ) 
-000E-00 3 
-000E-06 ) 


SHEEP = 


-O00E-06 ) 
.000E-06) 
-O00E-06 3} 
-000E-06) 
.000E-06 ) 
-O00E-06 } 
.O00E-06 } 


SWEEP = 


.O00E-06) 
.O00E-06 } 
.00CE-06 } 
-O00E-06 } 
-QO00E-06) 


8 4 


FOR 
FOR 
174% 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
175 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
176 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
177 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
FOR 
178 
FOR 
FOR 
FCR 
FOR 
FOR 


2.341E+10 
2.620E+10 


3.231E+04 
6.021E+06 
6.189E +06 
&.271E +06 
2.926E +08 
2.455E +10 
2.619E+10 


3.208E +04 
5.996E +06 
6.183E +06 
8.329E +06 
2.915E+08 
2.393E+¢10 
2.628E +10 


3.193E +04 
6.071E +06 
6.125£ +06 
8.259& +06 
2.909E +08 
2.455E +410 
2.628E+10 


3.235E&+04 
6.038E&+06 
6.128E +06 
8.291E +06 
2.934E+08 
2.466E +10 
2.622E 410 


3.272E +04 
6.095E& +06 
6.1156 +06 
&.309E +06 
3.010E+08 


oT 
TO0+* 3206 *2- 
TO+4S598°2- 
TO+ 3008" e- 
TO#3TS2Z°e- 
TO#3TL9°2- 
TO+*+32SS°2- 
TO#*3EL2°e- 
T0*3968°T- 
T0*4926°T- 
T0#3£90°T- 
00* 3662 °S- 
00+359T°9- 


$0-3082 °T- 
$0-362¢‘°T- 
$0-3089'T- 
$0-3220°T- 
90-3208°s- 
OT-3T80°T- 
OT-3589°T- 
OT-3086°9 
80-3598°T 


6 
20+3909'T- 
20*39STS°T- 
20*+3229°T- 
204#+3599 '‘[- 
20+3689°T- 
20#+3929°T- 
2043T29°T- 
20*+3920Z°T- 
2O+3b62ZL‘°T- 
2O0+3L22L°T- 
20*3022°T- 
20+3869'‘T- 


00+3622°- 
00+36TT‘S- 
00+3922°S- 
00+32S92°b- 
00+3229°s- 
00+322%° 2- 
00+306¢°T- 
(0=43Z2tS {- 
00+3T68°2 


8 
20+*+3S5S9¢ °T- 
20+329E°T- 
20+398°T- 
20+ 38TS°T- 
20#3095°T- 
20+3£09'°T- 
20*3969°T- 
20#+3929°T- 
2O#3AE2L°T- 
20*+399E2°T- 
20+362Z‘T- 
20#360Z°T- 


00+39828°S- 
00+ 3020°9- 
00+32ST *9- 
00+3059°S- 
00+ 3099 ° »- 
OO+AETE*<s- 
00+3968°T- 
T0-3950°2 

00+3899°¢ 


Z 
20#30S¢ ‘T- 
20+ 56S °T- 
2O+3TZ<e°T- 
20*+396E'°T- 
20+3ZTo‘T- 
20439065 '°T- 
20+350S'‘T- 
20#+3929S'°T- 
20+*+328S'T- 
20*328S'°T- 
20+3995S'‘T- 
2O+3ZUTS'T- 


00+3T8b* b- 
00+366S °5- 
00+3269° o- 
00+ 3582 °5- 
00+398S°s- 
00+39%6°2- 
00+3220°T- 
T0-3095°6 

00+3982°£ 


9 
20*+3292°T- 
20*+3892°T- 
20*3822°T- 
20+*3562‘T- 
20+*3905 °T- 
20+ 36TE‘*T- 
20+396¢ 'T- 
cO*3TZe°T- 
20+*3585 °T- 
20#49522°T- 
20¢39S6¢ °T- 
20+3TOs°T- 


OO+t3s62°T- TO-3TSb’s- Td-3S2Ee°Z 00+3255°T OO+3TZ2°T b =AT 
OO*3zsZs°T- T0-3806's- T0-3266°9 0O+3TT’T 00+3922'T S =AI 
00+*38Z5°T- T0-3206°b- [0-3229°S 00+36¢2°T 00+3650'T 9 =AI 
00+*3Z9=S°T- TO-322T°>- 1[0-3289°S OO+IeTe°T Oo+3ZTO'T ZL =AT 
00+36<0°T- TO-IJTOS°T- TO0-30To’*Z 00+3S5: °T 00+3590'°T @ =AI 
T0-3996°2- TO-3I22‘°d 00+36S2°T 00+3052‘°T 00+36S2°T 6 =AI 
T0-35952°2 OO+sece'T 00+*39%6°T OO+3TST°2 00+368¢°T OT =AI 
00+3620°2 00+3825°2 OOtFbEL°2 00+3899°2 00+390S5°T Tl =AI 
00+*39TS'Ss 00+3899'¢ 00+*]3098Z°£ 00+36¢2°¢ 00+3259°T 2T =AI 
TA JO S3INIWA 01314 

S z 2 T =xI 
COHZO9T*T- 2O*a2ZT'T- 2O+3aEGT°T- 20*39GT‘'IT- 20*302T°T- TT =AI 
ZCOCF69T° T= 2O*aedI*t- 2O+*aeet°T- 2O+JZ09T*I- 2O*392T°I- 2 =A 
COtZ6LT°T- Z2OHRGLT°I- 2O0*368T° T= 20+359T‘°T- 20*39S2T‘T- £ =AT 
COFCSTET’IT- 20*a98T°T- 2O*aZ6T°T- 20*¢3Z9T°T- 20¢322T°T- » =AT 
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OT+39Z9°2 OT+3L6E°2 GOFISTE*Z 904362E'9 90439E8T°9 9043966'S 4043902 20+3092'T 
OT+36T9°2 OT+3SSH°2 G049926'2 9O4FTL2°S 904368T'9 9043TZ20°9 HOFITE2°E 20+3002'T 
OT+3029°2 OT+3THS'2 GO43T9G*2 9043L0E°S 904356T'9 9043920°9 HOH35E2°= 20+30E2'T 
OT+39T9'2 OT+306E°2 G0439649°2 9043L0E'S 9043925T°9 9043650°9 40439902°E 20+3022°T 
OT+3T29°2 OL+3t6o°2 BO498S6°2 90497E2°S 904309T°9 9043066'F HOF306T'ES 20430TL'T 
OT+3TZ9°2 OT+36E6°2 BOHITIO'2 9OF9LTZ°S@ 9OF3LTI'9 9043500'9 04308T°E 20+3002°T 
OT+3L09°2 OT+39Sb°2 GO499TE'2 904392TZ2°S 9043SLT'9 9043896°S HOFI6TZ2°E 2043069'T 
OL+3909°2 OT+3STS°2 GO+3TH6'Z 904929T'S 9OF3EET'9 9043026°S HOF3TEZ*ES 2043089'T 
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APPENDIX Cc 


AIRFLOW CHARACTERIZATION SHEET 


TEST PROGRAM: Testing of Navy Aircraft in Arched Structure Hush House 


LOCATION: Kelly A.F.B 
DATE: _9 November 1982 


AIRCRAFT; F-4N 


ENGINE: ___J-79-CE-8 CELL DEPRESSION (in. wg.)s 1-7 
EVENT: , P) 


. CONDITIONS: Port Engine - A/B Power, Starboard Enpine = Idle Powe 
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AIRFLOW CHARACTERIZATION SHEET 


TEST PROGRAM: Thesis 

Location: _ Naval Air Station Jacksony 

DATE: NZA 

a ee : | 

ENGINE: J-79-GE-8 CELL DEPRESSION (ia. we.)s N/A 
N/A 


EVENT: —————————————————— 


Sinale Engine with A/B_@ cent rij 





CONDITIONS: 


‘AUGMENTER INLET 
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APPENDIX D 


EXHAUST FLOW DATA SHEET 


TEST PROGRAM; _7eSting of Navy Aircrafce in Arched Structure Hush House 


LOCATION: Kelly A.F.B. 9 San Antonio, Texas 


Oare; _9 November 1982 


AIRCRAFT; _ F-4-N 
ENGINE: J-79-CE-8 


EVENT: 2 
CONDITIONS: Port Engine - A/B, Starboard Engine = Idle Power 





AUGMENTER TUBE 


eterboard port 

















V=VELOCITY (feet/sec.) 


V= 209.5 V= 271.6 V= 239.7 
s ad @ 
T 454 T=454 





= 201.9 
. 


294.8 


AQ < 
tt 
4 @< 


454 454 


150.1 


454 


T=TEMPERATURE (degrees F) 


Vave= 2c 7ee 
Tavg= 454 


TOTAL MASS FLOW (tbm./eec.J)= 2067 


AUGMENTATION RATIO= 8.3 


eZ 


EXHAUST FLOW DATA SHEET 


TEST PROGRAM: _ Thesis 


Location: Naval Air Station, Jacksonville, Florida — 
DATE: N/A 


=) pz 
AIRCRAFT: & A 


Encine: J-/9-GE-8 
EVENT: N/A 


conoitions: cingle Engine with A/B @ centerline 


AUGMENTER TUBE 


eterboerd port 





V=VELOCITY (feet/sec.) 






ya 305.9 






s 
aA Gaye 







V= 318.3 
al 
T= 453.2 













v= 167.3 V=233 Mog 7 203 
ry e a 
T= 395.6 i= 331606 T= 395.6 






T=TEMPERATURE (degrees F) 


Vavg= 235.5 


FLOW (Ibm./eec.)= 2373 


AUGMENTATION RATIO“ 10.6 Tavg= 401.6 
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APPENDIX E 
DATA TRANSFER BY TAPE 

The IBM 3033 main-frame presently does not have "PHOTON" 
Graphics capabilities. In order to produce various visual 
plots that enhance the interpretation of tabulated data, the 
"DFOS" and "DF12" files must be properly formatted and then 
transferred to tape (for portability) allowing use of "PHOTON" 
on another computer which is “PHOTON" capable. 

The computer center will provide a non-labelled tape for 
Gata transfer. The operator, when queried by the programs 
below, will mount the tape allowing the "DFO9" and "DF12" 
files to be transferred to the tape. 

To prepare the files for transfer to tape a copy is made 
to ensure that the original data file is safe. Program "FORMO9 
JCL" properly formats and then writes the "DFO9" file to a new 
file called "DFO9C". Program "“FORM12 JCL" properly formats and 
then writes the "DF12" file to a new file called "“DF12C". The 
data sets "DFO9C" and "DF12C" which are produced using MVSHELP 
must be at least as large as the original "DFO9" and "DF12" 
data sets and created before running "FORM12 JCL" or “FORMO9 
JChN =. 

With the two files now properly formatted the file "DSTP 
JCL" is used to transfer the designated file to tape. In the 


following example the "DFO9C" file is transferred as the 


i, 3:0 


second file to the tape called "NICOLA". File "TSCAN JCL" then 
scans the tape to ensure the program was actually transferred 
to the tape. 

Two other files, "TAP JCL"” and "DUMP JCL", can dump a 
file from tape to the virtual machine and return the first 
block of the tape file respectively. With the files on tape 
and checked the tape can be dismounted by the operator at the 
computer center and transported to a "PHOTON" capable system 
in the "formatted" form. 


The various files referenced are now listed. 


FILE: DSTP JCL Al 


7/OSTP JOB (1541,9999), *ERIC',CLASS=E 

7/7 EXEC PGM=IEBGENER 

//7SYSPRINT OD SYSOUT=A 

“7SYSIN OD DUMMY 

//SYSUT1 DD DISP=SHR,DSN=MSS.S1541.PHOENICS.DFO9C 
//SYSUT2 DD UNIT =3400-6 ,VOL=SER=NICOLA ,DISP=(NEW,PASS), 


iA LABEL=(2,NL,,OUT), 
Cf DCB=( RECFM=FB ,LRECL=80 ,BLKSIZE=800 , DEN=3 ,OPTCD=Q) 
47 


Ike Dk 


FILE: TAP JCL Al 


//7TAP JOB (1541,9999),'TAP JCL',CLASS=E 

//7*MAIN CARDS=(50 ) 

// EXEC PGM=IEBGENER 

//SYSPRINT OD SYSOUT=A 

//7SYSIN DDO DUMMY 

//SYSUT1 0D UNIT=3400-6,VOL=SER=NICOLA ,DISP=(OLD,PASS ), 
17 LABEL=(2,NL,,IN), 


// OCB=(RECFM=FB,LRECL=80 ,BLKSIZE=800 ,OPTCD=Q ) 
//SYSUT2 00 SYSOUT=B,O0CB=BLKSIZE=80 
// 

FILE: DUMP JCL Al 


//7DUMP JOB (1541,9999), 'DUMP' ,CLASS=E 

// EXEC TAPE ,VOLIN=NICOLA 

//SYSIN OD * 
DMPREC(10,272,1),SKPFIL(10,1),XECUTE(999) 
/* 

4/ 


FILE: TSCAN JCL Al 


//GSCAN JOB (1541,9999), ‘NICOLA’ ,CLASS=E 
// EXEC TSCAN,VOLIN=NICOLA 
i 


2 


; 


FILE: FORMG9 JCL Al 


//FORMO9 JOB (1541,9999), CHAM’ ,CLASS=B 
/7*MAIN LINES=(20 ) 

// EXEC FORTVCLG 

/7FORT.SYSIN DD * 


gqganNnnga 


Oo 


an 


oe 


QANNANANNDNNANANNANAANANANDAN 


PROGRAM FORMO9 


Increase PARAMETER settings if more variables or a larger grid 
is required. 


PARAMETER( NPHIP=80 ) 

NB NXYMAX must be >= NZ and NX*NY 
PARAMETER( NXYMAX=20000 ) 
CHARACTER*4 NAMECNPHIP ) ,MESS( 10) 
CHARACTER*®15 FILIN,FILOUT 
CHARACTER*®1 ANS 
CHARACTER*®10 ACCESS 
LOGICAL LBUF( 10 ),STORE(NPHIP ) ,DEBUG 
DIMENSION RBUF(NXYMAX ),IBUF(10) 


100 FORMAT( 1X, 7110) 
101 FORMAT(6(1PE13.6)) 
102 FORMAT(1X%,79L1) 
103 FORMAT(1X%,1944 ) 


Reset DEBUG=.TRUE. to get error messages.. 
DEBUG=.TRUE. 


Set Logical Units... 
LUIN=9 
LUOUT=10 
LUBUG=6 


Set File Names... 

Name of unformatted restart file 
WRITE(6,%) * RESTART FILENAME?’ 
READ(5,‘(A)') FILIN 

Name of formatted restart file 
WRITE(6,%) ‘ FORMATTED FILENAME?‘ 
READ(5,‘'(A)') FILOUT 


Direct access or sequential? 

WRITE(6,%) ' DIRECT ACCESS (PHIDA) FILE? Y/N' 
READ(5,‘t(A1)') ANS 
IFC ANS.EQ.'Y'.OR.ANS.EQ.'Y*) THEN 

Version 1.4 default 
ACCESS='DIRECT' 
LENREC=256 
ELSE 

Version 1.3 
ACCESS='SEQUENTIAL' 


5, de 


FIte. 


OHOaNangaANDaADNH 


oO 


c 


FORMO9 JCL Al 
LENREC=100 
ENDIF 


WRITE(6,%) * DEFAULT LENREC? (100 V1.3, 256 V1.4) - Y/N‘ 
READ(5,°(A}*) ANS 

IFC ANS.NE.'Y*.AND.ANS.NE.‘°Y') THEN 

WRITE(6,%} ‘' LENREC? (WORDS NOT BYTES)‘ 

READ(5,‘'(I)*} LENREC 

ENDIF 


NBYTES is the number of bytes per single precision word. 


NBYTES=4 


Record length in bytes. 


LREC=NBY TES*¥LENREC 


C OPEN input file: 


C 


OPEN( UNIT=9, FORM=*UNFORMATTED' ,STATUS='‘OLD', 


1 ACCESS=ACCESS ,IOSTAT=IERR ,ERR=90 ) 


IFC IERR.NE.O) GO TO 90 
IF( ACCESS.EQ. ‘DIRECT') THEN 
READ( LUIN,REC=1)} (MESS(I),I=1,10) 
ELSE 
READ(LUIN) (MESS(I},I=1,10) 
ENDIF 
IF( DEBUG) WRITE( LUBUG,*) MESS 
IF(ACCESS.EQ. ‘DIRECT’ } THEN 
READ(LUIN,REC=2) (LBUF(I},I=1,4¢) 
ELSE 
READ(LUIN) (€LBUF(I},I=1,¢) 
ENDIF 
IF( DEBUG) WRITE( LUBUG,*) (LBUF(I),I=1,4¢) 
TFC ACCESS.EQ.*DIRECT’) THEN 
READ( LUIN,REC=3)} NX,NY »NZ»NPHI ,( IBUF(I),I=1,6),LENREC 
ELSE 
READ(LUIN) NX,NY>,NZ»NPHI,CIBUF(I },I=1,6),LENREC 
ENDIF 
NXNY =NX*¥NY 
IFENXNY .GT.NXYMAX}) GOTO 92 
IF(NPHI.GT.NPHIP} GO TO 93 
IF( DEBUG) WRITE( LUBUG,*) NX,NY NZ 


C OPEN output file: FORMATTED, SEQUENTIAL and of FIXED RECORD LENGTH. 


OPEN(UNIT=10,FORM='FORMATTED' ,STATUS='OLD', 

1 ACCESS='SEQUENTIAL ' ,IOSTAT=IERR,ERR=91 } 
IFCIERR.NE.O} GO TO 91 

WRITE( LUOUT,103) (MESS(I},I=1,10) 

WRITE(LUOUT,102) (LBUF(I},I=1,4)3 

WRITE( LUOUT,100)} NX,NY,NZ,NPHI,( IBUF(I),I=1,6),LENREC 

IFC ACCESS.EQ.‘*DIRECT'} THEN 

READ( LUIN,REC=4)} RINNER 

ELSE 

READ(LUIN}) RINNER 
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FILE: FORMO9 JCL Al 


ENDIF 
WRITE( LUOUT,101) RINNER 
IF( ACCESS.EQ. ‘DIRECT') THEN 
READ(LUIN,REC=5) (NAME(I ),I=1,NPHI ) 
ELSE 
READ(LUIN) (NAME(I ),I=1,NPHI ) 
ENDIF 
C Set C36... names for variables 51 to NPHI for version 1.4, to avoid 
C control characters etc] (Unless already stored) 
IFC ACCESS.EQ. ‘DIRECT’. AND.NPHI.GT.50) THEN 
READ(LUIN,REC=10) (STORE(I },I=1,NPHI ) 
DO 6 I=51,NPHI 
J=I-15 
IF(.NOT.STORE(I)) THEN 
IF(I.LT.100)} THEN 
WRITE(NAME(I),‘C''C'',I2,°° **)')V J 
ELSE 
WRITECNAME(I),‘'(''C’*,I3S)°) J 
ENDIF 
ENDIF 
6 CONTINUE 
ENDIF 
WRITE( LUOUT,103) (NAME( I )},I=1,NPHI ) 
IF( DEBUG) WRITE(LUBUG,*) (NAME(I),I=1,NPHI ) 
IF(ACCESS.EQ.‘DIRECT') THEN 
READ( LUIN,REC=6) (RBUF(I),I=1,NX) 
ELSE 
READ(LUIN) (RBUF(I },I=1,NX) 
ENDIF 
WRITE(LUOUT,101) (RBUF(I),I=1,NXx) 
IFCACCESS.EQ. ‘DIRECT') THEN 
READ(LUIN,REC=7)} (RBUF(I),I=1,NY } 
ELSE 
READ(LUIN) (RBUF(I},I=1,NY ) 
ENDIF 
WRITE( LUOUT,101) (RBUF(I),I=1,NY) 
IFC ACCESS.EQ. ‘DIRECT') THEN 
READ(LUIN,REC=8) (RBUF(I )},I=1,NZ) 
ELSE 
READ(LUIN) (RBUF(I ),I=1,NZ) 
ENDIF 
WRITE( LUOUT ,101) (RBUF(I)},I=1,N2Z) 
IF(ACCESS.EQ. ‘DIRECT’ ) THEN 
READ(LUIN,REC=9) (RBUF(I},I=1,NZ) 
ELSE 
READ(LUIN) (RBUF(I ),I=1,NZ) 
ENDIF 
WRITE( LUOUT,101) (RBUF(I),I=1,NZ) 
IF(ACCESS.EQ. ‘DIRECT'} THEN 
READ( LUIN,REC=10) (STORE(I),I=1,NPHI ) 
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FILE: FORMO9 JCL Al 


ELSE 
READ(LUIN) (STORE(1I),I=1,NPHI ) 
ENDIF 
WRITE€(LUOUT,102) (STORECI ),1=1,NPHI ) 
IF(DEBUG ) WRITE(LUBUG,102) (STORE(1I),ZI=1,NPHI ) 
NVAR=0 
DO 5 I=1,NPHI 
IFCSTORE(I)) NVAR=NVAR+1 
5 CONTINUE 
IREC=10 
DO 10 K=1,NZ 
DO 20 I=1,NVAR 
IFC ACCESS.EQ. *DIRECT') THEN 
IREC=IREC+1 
NVREC =NXNY/LENREC 
IF (MODO(NXNY »LENREC ).NE.0) NVREC=NVREC#+1 
DO 30 IVREC=1,NVREC 
JF=CIVREC-1 JXLENREC + 1 
JLSIVRECXLENREC 
JL=MINOC JUL, NXNY ) 
READ(LUIN,REC=IREC) (RBUF(J),J=JF,JL) 
30 CONTINUE 
ELSE 
READ(LUIN) (RBUF(J),J=1,NXNY ) 
ENDIF 
WRITE( LUOUT,101) (RBUF(J),J=1,NXNY ) 
20 CONTINUE 
10 CONTINUE 
IF( DEBUG ) WRITE(LUBUG,1900 ) 
1900 FORMAT(1X,' NORMAL STOP IN PROGRAM’ ) 


GO TO 999 
C Error trapping... 
90 IF( DEBUG) WRITE(LUBUG,1901) 
1901 FORMAT( 1X,‘ ERROR IN OPENING INPUT FILE‘) 
GO TO 999 
Ss] IF( DEBUG) WRITE(LUBUG,1911) 
1911 FORMAT(1X,‘ERROR IN OPENING OUTPUT FILE‘) 
GO TO 999 
Sz IFONXNY .GT.NXYMAX) WRITECLUBUG,1921) 


IFENZ.GT.NXYMAX ) WRITE( LUBUG,1923) 
1921 FORMAT(1X,'SETTING OF NX*®NY TOO LARGEs RESET IN FORMOS9 AND ' 
1,/7,1X,‘RE-CREATE THE LOAD MODULE ‘ ) 
1923 FORMAT(1X,'‘SETTING OF NZ TOO LARGES$ RESET IN FORMO9 AND ‘ 
1,/,1X,'RE-CREATE THE LOAD MODULE‘ ) 
GO TO 999 
93 WRITE ( LUBUG ,193 ) 
193 FORMAT(1X,‘SETTING OF NPHI TOO LARGE} RESET IN FORMOS AND '° 
1,/,1X,‘RE-CREATE THE LOAD MODULE‘) 
999 CONTINUE 
C CLOSE files ... 
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FILE: FORMO9 JCL Al 


CLOSE(LUIN,STATUS='‘KEEP ° } 

CLOSE( LUOUT ,STATUS='KEEP° ) 

STOP | 

END 
//GO.FTO9FOO1 DOD DSN=MSS.S1541.PHOENICS.DFO09,DISP=(OLD ) 
//GO.FT1OFOO1 DD DSN=MSS.S1541.PHOENICS.DFO9C ,DISP=(OLD ) 
/¥* 
// 


a7 


FILE: FORM12 JCL Al 


/7FORMI2 JOB (1541,9999), ‘CHAM’ ,CLASS=B 
/7*MAIN LINES=(20) 

// EXEC FORTVCLG 

//FORT.SYSIN DD * 


aNHaaggAaANg 


oO 


oO 


aAagaaaqanan 


File name: FORM12 last amended 15 October 1986 


PROGRAM FORM12 


cooce-- Program to rewrite a DF12 file into a formatted READCO file. 


(JPE & JFD ) 


Reset NXYIP if larger grid required.. 


PARAMETER( NXY1P=10000 ) 
DIMENSION AREC(NXYIP ) 


Reset fixed record length of DF12 file if required.... 


CHARACTER®30 FILIN,FILOUT 
LOGICAL DEBUG 

INTEGER SECTOR 

DATA LENREC/100/ 


100 FORMAT(315) 
101 FORMAT(5(1PE13.6)) 


SECTOR represents the largest record length that can be OPENed. 
SECTOR=65536 

NBYTES is the number of bytes per single precision word. 
NBYTES=4 
LREC=MINO(NBYTES*®LENREC ,SECTOR-NBYTES ) 
IBLK=LREC+NBYTES 


Set Logical Units... 
LUIN=12 
LUOUT=13 
LUBUG=6 

Reset DEBUG if error messages required.. 
DEBUG=.TRUE. 

Set file names 
FILIN='DF12 : 
FILOUT='DF12FORM ; 


OPEN input file: DIRECT ACCESS, UNFORMATTED and of FIXED 
RECORD LENGTH. 


OPEN( UNIT=12 ,FORM='UNFORMATTED' »STATUS=‘OLD', 
1 ACCESS='DIRECT',RECL=LREC, 
1 IOSTAT=IERR,ERR=90) 

IF(CIERR.NE.O) GO TO 90 


READ( LUIN,REC=1) NXP1,NYP1,NZP1 

NXY L=NXP1L¥NYP1 
IF(NXY1.GT.NXY1P) GO TO 92 

NREC=NXY1/LENREC 


JE ois) 


FILE: FORMI2 JCL Al 


ISHORT=NXY1-LENREC#NREC 
IFC ISHORT.GT.O) NREC=NREC+1 
IFCISHORT.EQ.0) ISHORT=LENREC 
IF(OEBUG) WRITE( LUBUG,*) NXP1,NYP1L,NZP1,NXY1,NREC ,ISHORT 
C 
C OPEN output file: FORMATTED, SEQUENTIAL and of FIXED RECORD LENGTH. 
C 
OPEN( UNIT=13 ,FORM='FORMATTED‘ ,STATUS=’OLD', 
1 ACCESS= ‘SEQUENTIAL ‘ ,IOSTAT=IERR,ERR=91 } 
IFCIERR.NE.O)} GO TO 91 
WRITE( LUOUT ,100 JNXP1,NYP1,NZP1 
NR=1 
DO § K=1,NZP1 
DO 6 IC=1,3 
IMAX=LENREC 
DO 7 JR=1,NREC 
IF(JUR.EQ.NREC } IMAX=ISHORT 
NR=NR+1] 
NRA=( JR-1 JXLENREC 
READ(LUIN,REC=NR) (AREC(NRA+I ),1=1,IMAX) 
7 CONTINUE 
c Write data to output file 
WRITE ( LUOUT ,101 JEAREC( II ),II=1,NXY1 ) 
6 CONTINUE 
5 CONTINUE 
IF(OEBUG) WRITE( LUBUG,1900 ) 
1900 FORMAT(1X, NORMAL STOP IN PROGRAM‘ } 
GO TO 999 
C Error trapping... 
90 IF( DEBUG ) WRITE( LUBUG,1901 ) 
1901 FORMAT(1X,'ERROR IN OPENING INPUT FILE‘) 
GO TO 999 
91 IF( DEBUG) WRITE( LUBUG,1911 3 
1911 FORMAT(1X,'ERROR IN OPENING OUTPUT FILE‘ ) 
GO TO 999 
92 IF( DEBUG) WRITE( LUBUG,1921) 
1921 FORMAT(1X, ‘GRID TOO LARGEs RESET NXY1P AND RE-CREATE', 
1 7,1X,'LOAD MODULE ' ) 
999 CONTINUE 
C CLOSE files ... 
CLOSE(LUIN,STATUS='KEEP' } 
CLOSE( LUOUT ,STATUS='KEEP' } 
STOP 
END 
//GO.FT12F001 DD ODSN=MSS.S1541.PHOENICS .DF12,DISP=( OLD) 
//GO.FT13F001 DD DSN=MSS.S1541.PHOENICS.DF12C ,DISP=(OLD } 
/* 
4/ 
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